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ABSTRACT 

We use cosmological, hydrodynamical simulations from the EAGLE and OWLS projects to 
assess the significance of recycled stellar ejecta as fuel for star formation. The fractional 
contributions of stellar mass loss to the cosmic star formation rate (SFR) and stellar mass 
densities increase with time, reaching 35% and 19%, respectively, at z = 0. The importance 
of recycling increases steeply with galaxy stellar mass for M* < 10 105 Mg, and decreases 
mildly at higher mass. This trend arises from the mass dependence of feedback associated 
with star formation and AGN, which preferentially suppresses star formation fuelled by recy¬ 
cling. Recycling is more important for satellites than centrals and its contribution decreases 
with galactocentric radius. The relative contribution of AGB stars increases with time and 
towards galaxy centers. This is a consequence of the more gradual release of AGB ejecta 
compared to that of massive stars, and the preferential removal of the latter by star formation- 
driven outflows and by lock up in stellar remnants. Recycling-fuelled star formation exhibits a 
tight, positive correlation with galaxy metallicity, with a secondary dependence on the relative 
abundance of alpha elements (which are predominantly synthesized in massive stars), that is 
insensitive to the subgrid models for feedback. Hence, our conclusions are directly relevant 
for the origin of the mass-metallicity relation and metallicity gradients. Applying the relation 
between recycling and metallicity to the observed mass-metallicity relation yields our best 
estimate of the mass-dependent contribution of recycling. For centrals with a mass similar to 
that of the Milky Way, we infer the contributions of recycled stellar ejecta to the SFR and 
stellar mass to be 35% and 20%, respectively. 
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formation 


1 INTRODUCTION 

The rate at which galaxies form stars is closely related to the 
amount of fuel that is available. Although we still lack a com¬ 
plete understanding of how galaxies obtain their gas, several po¬ 
tential sources of star formation fuel have been investigated in pre¬ 
vious works, both observationally and using hydrodynamical sim¬ 
ulations (e.g. Putman et al. 2009). Galaxies accrete gas from the 
intergalactic medium (IGM) along cold, dense, filamentary streams 
(e.g. Keres et al. 2005; Dekel et al. 2009; Brooks et al. 2009; van 
de Voort & Schaye 2012), which can extend far inside the halo 
virial radius, and through quasi-spherical infall from a diffuse hot 
halo, which contains gas that has been shock-heated to the halo 
virial temperature (Rees & Ostriker 1977; Silk 1977). Cosmolog¬ 
ical, hydrodynamical simulations give predictions for the relative 
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importance of these two ‘modes’ of gas accretion, generally indi¬ 
cating a dominant role for the cold mode in the global build up of 
galaxies, with the hot mode becoming increasingly important to¬ 
wards lower redshifts and in more massive systems (e.g. Birnboim 
& Dekel 2003; Keres et al. 2005, 2009; Crain et al. 2010; van de 
Voort et al. 2011; Nelson et al. 2013). Galaxies can also acquire 
new fuel for star formation by stripping the gas-rich envelopes of 
merging galaxies as soon as these become satellites in a group or 
cluster environment (e.g. Sancisi et al. 2008; van de Voort et al. 
201 1) or by re-accreting gas that has previously been ejected from 
the galaxy in an outflow and is raining back down in the form of a 
halo fountain (e.g. Oppenheimer & Dave 2008; Oppenheimer et al. 
2010). 

In addition to the various channels of accreting gas from the 
IGM, every galaxy has an internal channel for replenishing the 
reservoir of gas in the interstellar medium (ISM), namely the shed¬ 
ding of mass by the stellar populations themselves. Stars lose a 
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fraction of their mass in stellar winds before and while they go 
through the asymptotic giant branch (AGB) phase. Furthermore, a 
substantial amount of stellar material is released as stars end their 
lives in supernova (SN) explosions. Eventually, ~ 50% of the ini¬ 
tial mass of a stellar population will be released. If this material is 
not ejected into the circumgalactic medium (CGM), where it can 
emerge as X-ray emitting gas in the hot circumgalactic corona (e.g. 
Parriott & Bregman 2008; Crain et al. 2013), or entirely expelled 
from the galaxy into the IGM (e.g. Ciotti et al. 1991), but rather 
ends up in the cool ISM gas reservoir, then it can be ‘recycled’ to 
fuel subsequent generations of star formation (e.g. Mathews 1990; 
Martin et al. 2007). Note that what we call ‘gas recycling’ here 
is different from the process considered in works on galactic out¬ 
flow fountains, in which ‘recycling’ refers to the re-accretion of gas 
ejected from the ISM, regardless of whether it has ever been part 
of a star. In this work ‘recycled gas’ refers to the gas from evolved 
stars that is used to form new stars, regardless of whether it has 
been blown out of a galaxy. 

Using observational constraints on the rates of gas infall and 
the history of star formation, Leitner & Kravtsov (2011) assessed 
the significance of recycled stellar evolution products in the gas 
budget of a number of nearby disk galaxies (including the Milky 
Way). They modeled the global mass loss history of each galaxy 
from an empirically motivated distribution of stellar population 
ages and a set of stellar yields and lifetimes, and showed that the 
gas from stellar mass loss can provide most of the fuel required 
to sustain the current rates of star formation. They suggested that 
this internal supply of gas is important for fuelling star formation 
at late epochs, when the cosmological accretion rate drops or is 
suppressed by preventative feedback (e.g. Mo & Mao 2002), hence 
falling short of the observed star formation rate (SFR) of the galax¬ 
ies. Furthermore, Voit & Donahue (2011) argued that due to the 
high ambient pressures and the resulting short gas cooling times, 
central cluster galaxies are very efficient at recycling stellar ejecta 
into new stars. They showed that the stellar mass loss rates are gen¬ 
erally comparable to, or even higher than, the observed rates of star 
formation and emphasized the importance of including this form 
of internal gas supply in any assessment of the gas budget of such 
systems. These conclusions are consistent with the observation by 
Kennicutt et al. ( 1 994) that recycling of stellar ejecta can extend the 
lifetimes of gaseous discs by factors of 1.5 — 4, enabling them to 
sustain their ongoing SFRs for periods comparable to the Hubble 
time (see also Roberts 1963; Sandage 1986). These studies suggest 
that recycled stellar mass loss is an important part of the gas budget 
of star-forming galaxies, even hinting that it may be a necessary 
ingredient to reconcile the gas inflow and consumption rates of the 
Milky Way. 

In this paper, we investigate the importance of gas recycling 
for fuelling star formation by explicitly calculating the contribu¬ 
tion of stellar mass loss to the SFR and stellar mass of present-day 
galaxies. We use cosmological simulations from the Evolution and 
Assembly of GaLaxies and their Environments (EAGLE) project 
(Schaye et al. 2015, hereafter S15; Crain et al. 2015) to explore 
the recycling of stellar ejecta, as a cosmic average as a function 
of redshift and within individual (central and satellite) galaxies at 
z = 0, where we give quantitative predictions for recycling-fuelled 
star formation as a function of galaxy stellar mass and establish a 
connection with observational diagnostics by relating these predic¬ 
tions to gas-phase and stellar metallicities. 

The EAGLE simulations explicitly follow the mass released 
by stellar populations in the form of stellar winds and SN explo¬ 
sions of Types la and II, enabling us to study the relative signifi¬ 


cance of these mass loss channels for fuelling star formation. The 
subgrid parameters in the models for feedback associated with star 
formation and active galactic nuclei (AGN) have been calibrated to 
reproduce the 2 ~ 0 observed galaxy stellar mass function (GSMF) 
and the relation between stellar mass, M„, and the mass of the 
central supermassive black hole (BH), Mbh, with the additional 
constraint that the sizes of disc galaxies must be reasonable. The 
EAGLE simulations not only successfully reproduce these key ob¬ 
servational diagnostics with unprecedented accuracy, but are also 
in good agreement with a large and representative set of low- and 
high-redshift observables that were not considered during the cal¬ 
ibration (S15, Crain et al. 2015; Furlong et al. 2015; Lagos et al. 
2015; Rahmati et al. 2015; Sawala et al. 2015; Schaller et al. 2015a; 
Trayford et al. 2015). 

We consider the reproduction of a realistic galaxy population 
to be a prerequisite for this study, since its conclusions are sensitive 
to the detailed evolution of the gas ‘participating’ in galaxy forma¬ 
tion, requiring that the simulations accurately model the evolving 
balance between the inflow of gas onto galaxies and the combined 
sinks of star formation and ejective feedback. That EAGLE satis¬ 
fies this criterion is particularly advantageous, since hydrodynam- 
ical simulations are not subject to several limiting approximations 
inherent to simpler techniques, for example semi-analytic models 
of galaxy formation. This, in addition to their inclusion of a detailed 
implementation of chemodynamics, makes the EAGLE simulations 
an ideal tool for establishing quantitative predictions concerning 
the role of gas recycling in fuelling star formation. 

We also briefly explore the sensitivity of our results to the 
physical processes in the subgrid model. To do so, we use a suite of 
cosmological simulations from the Overwhelmingly Large Simu¬ 
lations (OWLS) project (Schaye et al. 2010). As the OWLS project 
aimed to explore the role of the different physical processes mod¬ 
elled in the simulations, it covers a wide range of subgrid imple¬ 
mentations and parameter values, including extreme variations of 
the feedback model and variations of the stellar initial mass func¬ 
tion (IMF). We will show that the efficiency of the feedback as¬ 
sociated with star formation and AGN plays an important role in 
regulating the fuelling of star formation with recycled stellar ejecta. 

We note that, because of the tight correlation we find between 
the contribution of stellar mass loss to the SFR (stellar mass) and 
the ISM (stellar) metallicity, our characterization and explanation 
of the role of stellar mass loss as a function of galaxy mass and type 
has important and direct implications for the origin of the mass- 
metallicity relation. 

This paper is organized as follows. In Section 2 we present 
a brief overview of the simulation set-up and the subgrid modules 
implemented in EAGLE. In this section we also introduce the two 
quantities we use to assess the importance of gas recycling, namely 
the fractional contributions of stellar mass loss to the SFR and stel¬ 
lar mass. In Section 3 we present quantitative predictions from EA¬ 
GLE for the evolution of the cosmic averages of these quantities 
and for their dependence on metallicity and galaxy stellar mass. We 
explore the sensitivity of these results using a set of OWLS simula¬ 
tions in Section 4. Finally, we summarize our findings in Section 5. 


2 SIMULATIONS 

The amount of gas that galaxies can recycle to form new genera¬ 
tions of stars, depends fundamentally on the fraction of stellar mass 
that is returned to the ISM. How much of this mass is actually used 
to fuel star formation is not straightforward to calculate analyti- 
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cally, due to the variety of processes, such as cosmological infall, 
gas stripping of satellite galaxies, and feedback associated with star 
formation and AGN, that can have an effect on the star formation 
histories of individual galaxies. Hence, we use cosmological simu¬ 
lations from the EAGLE and OWLS projects to investigate this. 

For the majority of this work we use the EAGLE sim¬ 
ulations, which were run with a modified version of the 
smoothed particle hydrodynamics (SPH) code GADGET3 (last 
described by Springel 2005) using a pressure-entropy formu¬ 
lation of SPH (Hopkins 2013; see Schaller et al. 2015b for 
a comparison between SPH flavours). The simulations adopt a 
ACDM cosmology with parameters [fi m , fib, Ha, 0's, n s , h] = 
[0.307, 0.04825, 0.693, 0.8288,0.9611, 0.6777] (Planck Collabo¬ 
ration et al. 2014). 

We will study primarily the largest EAGLE simulation, which 
we will refer to as Ref-L100N1504 (as in S15) or as the ‘fidu¬ 
cial’ model. This simulation was run in a periodic volume of size 
L = 100 comoving Mpc (cMpc), containing N = 1504 3 dark 
matter particles and an equal number of baryonic particles. The 
gravitational softening length of these particles is 2.66 comoving 
kpc (ckpc), limited to a maximum physical scale of 0.7 proper kpc 
(pkpc). The particle masses for baryons and dark matter are initially 
m b = 1.8 x 10 6 Mq and mjm = 9.7 x 10 6 M@, respectively. 
However, during the course of the simulation the baryonic particle 
masses change as mass is transferred from star to gas particles, cor¬ 
responding to the recycling of mass from stellar populations back 
into the gas reservoir. 


2.1 Subgrid physics 

The subgrid physics used in EAGLE is largely based on the set of 
subgrid recipes developed for OWLS, but includes a few impor¬ 
tant improvements. Star formation is modelled using a metallicity- 
dependent density threshold (given by Schaye 2004), above which 
gas particles are assigned a pressure-dependent SFR (that by con¬ 
struction reproduces the observed Kennicutt-Schmidt star forma¬ 
tion law; Schaye & Dalla Vecchia 2008) and are converted stochas¬ 
tically into star particles. Each star particle represents a stellar pop¬ 
ulation of a single age (simple stellar population; SSP) and inher¬ 
its its mass and metallicity from its progenitor gas particle. The 
adopted IMF is a Chabrier (2003) IMF, spanning the mass range 
of 0.1 — 100 Mq. Following the prescriptions of Wiersma et al. 
(2009b), an SSP loses mass through stellar winds and supernova ex¬ 
plosions (SN Type II) from massive stars and through AGB winds 
and SN Type la explosions from intermediate-mass stars. The time- 
dependent mass loss, which we show in Section 2.2, is calculated 
using the metallicity-dependent stellar lifetime tables of Portinari 
et al. (1998), in combination with the set of nucleosynthetic yields 
of Marigo (2001) (for stars in the mass range 0.8 — 6 Mg ) and 
Portinari et al. (1998) (for stars in the mass range 6 — 100 Mq), 
all of which are based on the same Padova evolutionary tracks. For 
SN Type la, the yields are taken from the W7 model of Thiele¬ 
mann et al. (2003) and the distribution of progenitor lifetimes is 
modelled using an empirically motivated time-delay function that 
is calibrated to reproduce the observed cosmic SN Type la rate 
(see fig. 3 of S15). At every gravitational time step (every 10th 
time step for star particles older than 100 Myr), the ejecta are dis¬ 
tributed over the neighbouring gas particles according to the SPH 


interpolation scheme 1 . The simulations follow the abundances of 
11 individual elements, which are used to calculate the rates of ra¬ 
diative cooling and heating on an element-by-element basis and in 
the presence of Haardt & Madau (2001 ) UV and X-ray background 
radiation (Wiersma et al. 2009a). Energy feedback from star for¬ 
mation is implemented by stochastically injecting thermal energy 
into the gas surrounding newly-formed star particles as described 
by Dalla Vecchia & Schaye (2012). The fraction f t h of the total 
available feedback energy that is used to heat the gas, depends on 
the local gas metallicity and density, so as to account for increased 
thermal losses in higher metallicity gas and to compensate for the 
increased numerical radiative losses in higher density gas (Crain 
et al. 2015). The growth of BHs is modelled by inserting seed BHs 
into haloes more massive than mhaio.mm = 10 10 ft -1 M 0 , which 
can grow either through gas accretion, at a rate that depends on the 
angular momentum of the gas, or through mergers with other BHs 
(Booth & Schaye 2009; Rosas-Guevara et al. 2015). AGN feedback 
is implemented as the stochastic injection of thermal energy into the 
gas surrounding the BH (Booth & Schaye 2009; Dalla Vecchia & 
Schaye 2012). The subgrid routines for stellar and AGN feedback 
have been calibrated to reproduce observations of the present-day 
GSMF. the Al* — Mbh relation and to yield reasonable galaxy sizes 
(S15, Crain et al. 2015). 

2.2 Mass released by an SSP 

Fig. 1 shows the total mass (left panel) and metal mass (right 
panel) released by an SSP as a function of its age as prescribed by 
the chemodynamics model. The curves show the total integrated 
mass ejected (black) and the integrated mass sourced by AGB stars 
(blue), massive stars (i.e. stellar winds plus SN Type II; purple) 
and SN Type la (cyan) for two different SSP metallicities: solar 
(solid lines) and 1 percent of solar (dashed), using a solar value 
of Z 0 = 0 .0127. Both panels show that the ejected (metal) mass, 
which is expressed as a fraction of the total initial mass of the SSP. 
increases as the SSP ages. Initially only massive stars contribute, 
but for ages > 10 8 yr the contribution from AGB stars becomes 
increasingly significant. Comparing, for each channel, the total to 
the metal mass loss shows that the ejecta from massive stars are 
more metal-rich than those from AGB stars. The contribution from 
SN Type la to the (metal) mass loss remains insignificant for all 
SSP ages' 2 . Varying the metallicity over two orders of magnitude 
changes the total mass loss by only a few percent, but changes the 
total ejected metal mass, as well as the relative contribution from 
massive stars, by ~ 10 — 15%. 

Since the choice of IMF determines the relative mass in 
intermediate-mass and massive stars per unit stellar mass formed, 
it affects the mass loss from an SSP. Leitner & Kravtsov (2011) 
indeed show that the differences between alternative, reasonable 

1 As discussed in S15 and different from Wiersma et al. (2009b), EAGLE 
uses weights that are independent of the current gas particle mass for the 
distribution of stellar mass loss. The reason for this is to avoid a runaway 
process, causing a small fraction of the particles to end up with very large 
masses compared to their neighbours, as particles that have grown massive 
due to enrichment, are also likely to become increasingly enriched in future 
time steps, if they cany more weight in the interpolation. 

2 Note that we only show the relative contributions from massive and 
intermediate-mass stars to the total ejected metal mass. These may be dif¬ 
ferent from their contributions to the ejected mass of individual elements, as 
for example iron, which has a substantial fraction of its abundance sourced 
by SN Type la (see fig.2 of Wiersma et al. 2009b). 
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Figure 1 . The cumulative fraction of the initial mass (total: left panel; in the form of metals, i.e. elements heavier than helium: right panel) that is released by 
an SSP as a function of its age, adopting a Chabrier (2003) IMF in the range 0.1 — 100 Mg. The curves show the contributions from AGB stars (blue), massive 
stars (purple) and SN Type la (cyan), as well as the total (metal) mass ejected by the SSP (black), for two stellar metallicities: solar (solid) and 1 percent of 
solar (dashed). Initially, only massive stars contribute to the mass loss, but for SSP ages > 10 s yr the contribution from AGB stars becomes increasingly 
significant. These AGB ejecta are, however, less metal-rich than the ejecta from massive stars. The contribution from SN Type la to the (metal) mass loss 
remains insignificant at all times. Increasing the metallicity does not have a strong effect on the total mass loss, but increases the total ejected metal mass as 
well as the relative contribution from massive stars. 


choices of the IMF can be significant (see their fig. 1). In Ap¬ 
pendix A, we similarly conclude that the total and metal mass loss 
is a factor of ~ 1.5 greater for a Chabrier IMF than for the more 
bottom-heavy Salpeter IMF (which is adopted by one of the OWLS 
model variations examined in Section 4). 


2.3 Numerical convergence 

In order to test for numerical convergence, we use a set of three 
simulations that were run in volumes of size L = 25 cMpc. This 
includes a high-resolution simulation ( Recal-L025N0752 ), whose 
subgrid feedback parameters were recalibrated to improve the fit 
to the observed present-day GSMF (see table 3 of SI5). We show 
a concise comparison between the fiducial simulation and Recal- 
L025N0752 when we present results as a function of halo and stel¬ 
lar mass in Section 3.3, while a more detailed convergence test can 
be found in Appendix Bl. In the rest of the results section (Sec¬ 
tion 3) we use only the fiducial simulation, which, due to its 64 
times greater volume than Recal-L025N0752 , provides a better sta¬ 
tistical sample of the massive galaxy population, and models a more 
representative cosmic volume. 

2.4 Identifying haloes and galaxies 

Haloes are identified using a Friends-of-Friends (FoF) algorithm 
(Davis et al. 1985), linking dark matter particles that are separated 
by less than 0.2 times the mean interparticle separation. Gas and 
star particles are assigned to the same halo group as their nearest 
dark matter particle. The SUBFIND algorithm (Springel et al. 2001; 
Dolag et al. 2009) then searches for gravitationally bound substruc¬ 


tures within the FoF haloes, which we label ‘galaxies’ if they con¬ 
tain stars. The galaxy position is defined to be the location of the 
particle with the minimum gravitational potential within the sub¬ 
halo. The galaxy at the absolute minimum potential in the FoF halo 
(which is almost always the most massive galaxy) is classified as 
the ‘central’ galaxy, whereas the remaining subhaloes are classified 
as ‘satellite’ galaxies. 


The mass of the main halo, M 200 , is defined as the mass in¬ 
ternal to a spherical shell centred on the minimum gravitational po¬ 
tential, within which the mean density equals 200 times the critical 
density of the Universe. The subhalo mass, M su b, corresponds to 
all the mass bound to the substructure identified by SUBFIND. The 
stellar mass, M*, refers to the total mass in stars that is bound to 
this substructure and that resides within a 3D spherical aperture of 
radius 30 pkpc. Other galaxy properties, such as the SFR and the 
stellar half-mass radius, are also computed considering only parti¬ 
cles within this aperture, mimicking observational measurements of 
these quantities (as shown in fig 6. of S15, the present-day GSMF 
using a 30 pkpc 3D aperture is nearly identical to the one using 
the 2D Petrosian aperture applied by SDSS). The aperture has neg¬ 
ligible effects on stellar masses for M* < 10 11 Mq and galactic 
SFRs, as the vast majority of the star formation takes place within 
the central 30 pkpc. For the more massive galaxies, on the other 
hand, the stellar masses are somewhat reduced, as the aperture cuts 
out the diffuse stellar mass at large radii that would contribute to 
the intracluster light. 
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2.5 Measuring the star formation rate and stellar mass 
contributed by recycling 


We explicitly track the contributions to the SFR and stellar mass 
from gas recycling. For a gas particle of mass rn g (t) at time t, the 
total fraction of its mass contributed by released stellar material (in 
the form of hydrogen, helium and heavy elements) is given by 


/g,rec(f) = 


m’(t) — mb 


(1) 


where mb is the initial gas mass of gas particles at the start of the 
simulation. Since a gas particle is the smallest quantum of mass we 
are able to consider, its recycled fraction is by construction assumed 
to be perfectly mixed. Therefore, if the gas particle is considered 
star-forming, /g, rec (f) also indicates the fraction of its current SFR 
that is contributed by stellar ejecta. Then, summing up the contribu¬ 
tions from all iVf al gas particles in a galaxy (within the 30 pkpc 3D 
aperture) yields the SFR contributed by recycling for this galaxy: 


AfJ al 


SFR 


gal 


w = E 


mj(i) — mb 
m|(t) 


SFR‘(t), 


( 2 ) 


where SFR l (t) is the SFR of gas particle i at time t. Similarly, 
summing up the contributions from all N g ° B gas particles in the 
simulation volume yields the cosmic average of this quantity: 


SFR^it) = Y mb SFR j t y 


(3) 


Since a star particle inherits its mass and elemental abun¬ 
dances from its progenitor gas particle, the fraction of its mass con¬ 
tributed by recycling is: 




K.init - m b 


*,init 


(4) 


where ?n 4 init = m|(fbirth) is the mass of star particle j at the 
time of its birth, fbirth- Note that equation (4) is valid for all t ^ 
t birth, even though the star particle itself loses mass. This is again 
a consequence of the assumption of perfect mixing on the particle 
scale. Summing up the contributions from all iVf al star particles in 
a galaxy that are within the 3D aperture, 


7 V gal 


Ms:Ut) = E 


;,i„it - 


!(*). 


j=l 

and all N^ OB star particles in the simulation volume, 


M vL, mi ■ — mb 

cos /,\ \^ *,imt u 

*,rec(^7 / . j ^ 


m, 


j=i 


*,init 


(5) 


( 6 ) 


give the galaxy stellar mass and cosmic stellar mass, respectively, 
contributed by recycled gas. 

While SFR lec and M*, re c are related, it is still helpful to con¬ 
sider both: SFR Tec indicates the instantaneous impact of gas recy¬ 
cling, whereas M„, re c indicates the importance of recycling over 
the past history of star formation. In this work we mainly focus 
on the relative contribution of gas recycled from stellar mass loss 
to the total (cosmic or galactic) SFR and stellar mass. Normaliz¬ 
ing SFR lec and A/*, rec by the respective total quantities, yields 
SFR lec /SFR and M* !rec /M*, specifying the fractions of the SFR 
and the stellar mass that are due to stellar mass loss. 

In addition to the total amount of recycling, we will also con¬ 
sider the relative contributions from the different sources of stellar 


mass loss that were included in the subgrid model (Section 2.2). As 
the transfer of mass from AGB stars, SN Type la and massive stars 
between star and gas particles is explicitly followed by the EA¬ 
GLE simulations 5 , we can calculate SFR lec /SFR and M*, re c/M» 
solely due to gas from AGB stars by simply replacing m g (t) — mb 
in equations (2) and (3) by m\ GB and replacing m? t init — mb in 
equations (5) and (6) by m J AGB init , where ttlagb is the mass from 
AGB stars in the respective gas or star particle. The SFR lec /SFR. 
and Mt,iec/M t due to gas from SN Type la and massive stars are 
calculated analogously. 


3 RECYCLED STELLAR MASS LOSS IN EAGLE 

In this section we use the fiducial EAGLE simulation, Ref- 
L100N1504, to make quantitative predictions for the importance 
of gas recycling for fuelling ongoing star formation in present-day 
galaxies over a wide range of galaxy masses. However, we start 
with a brief investigation of the evolution of recycling-fuelled star 
formation over cosmic history. 

3.1 Evolution of the cosmic average 

The left panel of Fig. 2 shows the total cosmic SFR density (black), 
the cosmic SFR density fuelled by stellar mass loss (red, solid: ‘re¬ 
cycled’) and the cosmic SFR density fuelled by unprocessed gas 
(blue: ‘non-recycled’) as a function of redshift. The red curve has 
been split into the contributions from the three mass loss channels 
that are tracked by the simulation: AGB stars (dashed), massive 
stars (dotted) and SN Type la (dot-dashed). To get a better idea of 
the evolution of th & fractional contribution from recycled gas to the 
cosmic star formation history, we show the evolution of the cos¬ 
mic average SFR rec /SFR, as well as the fractional contribution 
per channel, in Fig. 3 (red). 

At z > 2 there is little difference between the SFR density 
due to ‘non-recycled’ gas and the total SFR density. At these high 
redshifts most of the fuel for star formation is due to unprocessed 
gas 4 , since there has simply not been much time for stellar popula¬ 
tions to evolve and to distribute a significant amount of gas that can 
be recycled. From the ‘recycled’ curve we see that the SFR density 
fuelled by recycled stellar mass loss rises rapidly at high redshift, 
peaks at z ~ 1.3, and then declines steadily towards 3 = 0. This 
trend is similar to the evolution of the total SFR density, although 
with a delay of ~ 1.5 Gyr (the total SFR density peaks at z ~ 2). 
Furthermore, the slope of the ‘recycled’ curve is steeper at high red¬ 
shift and shallower at low redshift compared to that of the total SFR 
density, indicating that gas recycling becomes increasingly impor¬ 
tant for fuelling star formation. This is consistent with the rapid 
rise of the total SFR IOC /SFR with decreasing redshift in Fig. 3. 
Our fiducial EAGLE model indicates that 35% of the present-day 
cosmic SFR density is fuelled by recycled stellar mass loss. 

The right panel of Fig. 2 shows the build up of the cosmic 
stellar mass density, the total as well as the contributions from re¬ 
cycled and unprocessed gas. The evolution of the cosmic average 

3 Note that these enrichment channels only refer to the last enrichment 
episode. Every stellar population releases mass via the different channels in 
a way that depends only on its age and metallicity (for a given IMF). 

4 Note that this does not imply that most of the SFR, and hence stellar 
mass, is in the form of Pop III (i.e. metal-free) stars, because the stellar evo¬ 
lution products are mixed with the unprocessed material (in the simulations 
on the scale of a gas particle). 
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Figure 2. The evolution of the cosmic SFR density (left) and the cosmic stellar mass density (right) fuelled by recycled stellar mass loss (red), as well as 
the SFR and stellar mass densities fuelled by all gas (black) and gas that has not been recycled (blue). The ‘recycled’ SFR and stellar mass densities are 
split according to the contributions from AGB stars (dashed), massive stars (dotted) and SN Type la (dot-dashed). Recycling of stellar mass loss becomes 
increasingly important for fuelling star formation towards the present day. The gas from massive stars accounts for the majority of the cosmic SFR and stellar 
mass density from recycled gas at high redshift, but the contribution from AGB stars increases with time (accounting for the majority of the ‘recycled’ SFR 
density for z < 0.4). 


M*, rec /M* is shown in Fig. 3 (blue). The stellar mass density is 
related to the SFR density, as one can calculate the former by inte¬ 
grating the latter over time (while taking into account stellar mass 
loss). Hence, similar to the SFR density, the stellar mass density 
is initially (z > 2) dominated by star formation from unprocessed 
gas, while the contribution from recycling becomes increasingly 
important towards z = 0. EAGLE indicates that, at the present 
day, 19% of the cosmic stellar mass density has been formed from 
recycled stellar mass loss. 

Comparing the different sources of stellar mass loss, we see 
that massive stars initially account for the majority of the SFR and 
stellar mass density from recycled gas. These stars have short life¬ 
times and are therefore the first to contribute to the mass loss from a 
stellar population (see Fig. 1). Towards lower redshift the mass lost 
by AGB stars becomes increasingly important and even becomes 
the dominant contributor to the SFR density from recycled gas for 
z < 0.4 (while remaining subdominant in the case of the stellar 
mass density). As expected from Fig. 1, recycled SN Type la ejecta 
do not contribute significantly to the cosmic SFR density at any 
redshift. 

3.2 Relation with metallicity at z = 0 

Having studied the evolution of the cosmic average SFR Tec /SFR 
and M*, rec /AF*, we will now take a closer look at the z = 0 val¬ 
ues for individual galaxies in the Ref-L100N1504 simulation. In 
the next section we will give predictions for the fuelling of star for¬ 
mation by recycled stellar ejecta in present-day central and satellite 
galaxies as a function of their halo and stellar mass. To be able to re¬ 
late these predictions to observational diagnostics, we first explore 
the relation between recycling-fuelled star formation and present- 
day metallicity. We will show that the fact that metals are synthe- 



Figure 3. The evolution of the fractional contribution of recycled stellar 
mass loss to the cosmic SFR density (red) and cosmic stellar mass density 
(blue), where we show the total (solid) as well as the contributions from 
AGB stars (dashed), massive stars (dotted) and SN Type la (dot-dashed). 
With decreasing redshift, an increasing fraction of the cosmic SFR and stel¬ 
lar mass density is fuelled by recycled gas, which we find to be 35% and 
19%, respectively, at z = 0. 
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Figure 4. The fractional contribution of recycled stellar mass loss to the SFR (left) and stellar mass (right) of central galaxies at z = 0 as a function of their 
average ISM and stellar metallicity, respectively. The grey scale indicates the number of galaxies in each cell, where we only include galaxies with stellar 
masses corresponding to at least 100 gas particles. In the left panel we only consider subhaloes with a non-zero SFR. We find tight power-law relations 
between the recycled gas contributions and the respective metallicity measures. These relations exhibit a slight mass dependence as a result of the increasing 
contribution from massive stars relative to intermediate-mass stars to the SFR and stellar mass for M* > 10 ll) o Mg. The best-fit relations (equations 10 and 
1 1), plotted for galaxies with M, ~ 10 9 ' 5 Mq (red, solid line), M* ~ 10 10 ' 5 Mg (blue, dot-dashed line) and M* ~ 10 11 ' 5 Mg (purple, dashed line), 
enable one to estimate the importance of gas recycling in present-day galaxies from their observed metallicity and cr-enhancement. 


sized in stars and are distributed over the ISM as the evolving stel¬ 
lar populations lose mass, makes them an excellent observational 
proxy for the contribution of stellar ejecta to the SFR and stellar 
mass. 

To study the SFR lec /SFR, we only consider subhaloes with 
a non-zero 5 SFR, while to study the A/», rec /Af*, we only consider 
subhaloes with a non-zero stellar mass. For our fiducial simulation 
this yields samples of 44 248 and 325 561 subhaloes, respectively. 
In this section, however, we additionally require the subhaloes to 
have a galaxy stellar mass corresponding to at least 100 gas parti¬ 
cles, which yields samples of 14 028 and 16 681 subhaloes, respec¬ 
tively. 

Fig. 4 shows the fraction of the SFR (left panel) and stellar 
mass (right panel) fuelled by recycling as a function of, respec¬ 
tively, the mass-weighted absolute metallicity Z gas of ISM gas (i.e. 
star-forming gas) and the mass-weighted absolute metallicity Z* 
of stars, both for present-day central galaxies 6 * . We find strong cor¬ 
relations between these quantities, with more metal-rich galaxies 
having a larger fraction of their SFR and stellar mass contributed 
by recycling. The figure reveals tight power-law relations between 
SFR lec /SFR and Z gas , characterized by a Pearson correlation co¬ 
efficient of 0.95, and between M„ jrec /M* and Z„, with a cor¬ 
relation coefficient of 0.99. For the former we find a lcr scatter 
of ~ 0.1 — 0.2 dex for Z gas < 10 -1 ' 9 and < 0.05 dex for 
Z gas > 10 -1 ' 9 , while for the latter we find an even smaller lcr scat¬ 
ter of ~ 0.01 — 0.03 dex. Furthermore, as we show in Appendix B2, 

5 ‘Non-zero’ means containing at least one star-forming gas particle, which 
corresponds to a specific SFR (= SFR/M *) of > 10 -12 yr -1 at M* ~ 
10 9 Mg and > 10“ 14 yr” 1 at M* ~ 10 11 Mg. 

6 Although we do not explicitly show it, the results for central galaxies 

presented in this section are consistent with the results for satellite galaxies. 


both relations are converged with respect to the numerical resolu¬ 
tion. 

The tight relation between the contribution of recycled gas 
to star formation and metallicity is not surprising considering that 
heavy elements were produced in stars and that their abundance 
must therefore correlate with the importance of stellar ejecta as star 
formation fuel. The contribution of recycling to the stellar mass is 

equal to the ratio of the mean stellar metallicity (^ Z* ^) and the 

mean metallicity of the ejecta ( Z e j ^) that were incorporated into 
the stars, 



The same holds for the contribution of stellar mass loss to the SFR. 



The metallicity of the ejecta depends on the age and metallicity 
of the SSP, as well as on the IMF. From Fig. 1 we can see that 

for our (Chabrier) IMF, (^Z e ^j « 0.033/0.45 « 0.073 for a 10 

Gyr old SSP with solar metallicity. Hence, log 10 M* irec /M* « 
log 10 Z* + 1.1, where the slope and normalization are close to the 
best-fit values that we determine below. Note that using ages of 100 
Myr and 10 Myr instead of 10 Gyr gives normalizations of 0.91 and 
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Figure 5. The a-element-to-iron abundance ratio of central galaxies at 
z = 0 as a function of stellar mass. We show [a/Fe], represented by 
[O/Fe] (solid) and [Mg/Fe] (dotted), of ISM gas (green) and stars (blue) 
as predicted by EAGLE, and compare with observations of the stellar 
[a/Fe] from Thomas et al. (2010) (converted to a solar abundance ratio of 
X ( - > /A'g e = 4.44). The curves show the median value in each logarithmic 
mass bin of size 0.2 dex, if it contains at least 10 haloes and the stellar mass 
corresponds to at least 100 gas particles. The shaded regions mark the 10th 
to 90th percentiles, shown only for [O/Fe]. For M* < 10 10,5 Mg, [O/Fe] 
([Mg/Fe]) is approximately constant at ~ 0.1 (—0.2) for gas and at 0.25 
(—0.05) for stars. For M* > 10 10 5 Mg, [O/Fe] and [Mg/Fe] increase 
with stellar mass, in such a way that the slope matches the observations, 
reflecting the enhancement in the contribution to the SFR and stellar mass 
from massive stars relative to that from intermediate-mass stars. 


0.77, respectively. Using an age of 10 Gyr but a stellar metallicity 
of 0.01 Zq instead of Zg yields a normalization of 1.2. 

There is, however, an additional factor at play that may dis¬ 
tort the one-to-one correlation between the contribution of recycled 
gas to the SFR (and therefore to the stellar mass) and metallicity, 
namely the relative significance of the different mass loss channels. 
This depends on the timescale on which stars are formed, but is also 
affected by processes like stellar and AGN feedback. Given that the 
ejecta from massive stars have ~ 4—6 times higher metallicity than 
those from intermediate-mass stars (dependent on metallicity; see 
Fig. 1), a higher contribution of the mass loss from massive stars to 
the SFR (for fixed SFR lec /SFR) would yield a higher ISM metal¬ 
licity, and would hence change the relation between SFR lec /SFR 
and Z gas . As we will show in Section 3.3.3, the contribution to the 
SFR of the mass loss from massive stars relative to that from AGB 
stars varies as a function of stellar mass, and in particular increases 
at the high-mass end. This introduces a mild mass dependence in 
the SFR rec /SFR- Z gas and M*, rec /M,- Z * relations 7 . In order 

' Another factor is that the metal yields depend on metallicity (Fig. 1). 
This can change the SFR Iec /SFR- Z gas and M*, rec /M«- Z* relations 
even if the contributions from the different channels remain fixed. However, 
even a factor of 100 variation in the metallicity changes the metallicity of 
the stellar ejecta by only a few percent, which is significantly smaller than 
the effect of the change in the relative channel contributions in massive 
galaxies. 


to relate this variation of the relative contribution from different 
mass loss channels to an observational diagnostic, we consider the 
average a-enhancement, [a/Fe], represented by [O/Fe] (as oxy¬ 
gen dominates the a-elements in terms of mass fraction), of ISM 
gas and stars. The fact that a-elements are predominantly synthe¬ 
sized in massive stars, whereas of iron ~ 50% is contributed by 
intermediate-mass stars in the form of SN Type la explosions and 
winds from AGB stars (e.g. Wiersma et al. 2009b), makes [a/Fe] 
a good tracer for the relative importance of massive stars. 

Adopting the usual definition of the abundance ratio, 

= logic - logic » (9) 

where X x is the mass fraction of element x and Xg/A'g e = 4.44 
is the solar abundance ratio (Asplund et al. 2009), we show [O/Fe] 
as a function of stellar mass in Fig. 5. The curves show the me¬ 
dian in logarithmic mass bins of size 0.2 dex that contain at least 
10 haloes and correspond to a stellar mass of at least 100 gas parti¬ 
cles. The shaded regions mark the 10th to 90th percentile ranges. In 
both the gas-phase (green, solid) and the stellar phase (blue, solid), 
[O/Fe] is approximately constant at ~ 0.1 and ~ 0.25, respec¬ 
tively, for M* < 10 10,5 Mg, but increases with stellar mass for 
M* > 10 10 ' 5 Mg. Comparing this to observations of the stellar 
[a/Fe] for a sample of 3360 early-type galaxies from Thomas et al. 
(2010) (best-fit relation, after correcting for the difference in the 
set of solar abundances used; black dashed line), we find excellent 
agreement in terms of the slope and the normalization. While this is 
encouraging, suggesting that we capture the right mass dependence 
in the SFR vec /SFR- Z gas and M*, rec /M«- Z* relations and that 
the cooling rates (which are dominated by oxygen at T ~ 2 x 10° 
K and by iron at T ~ 10° K; see Wiersma et al. 2009a) employed 
by the simulation are realistic, the predicted abundance ratio is un¬ 
certain by a factor of > 2 due to uncertainties in the nucleosyn- 
thetic yields and SN Type la rate (Wiersma et al. 2009b). It is there¬ 
fore somewhat surprising that the agreement in the normalization 
is this good. If we consider [Mg/Fe], which is another indicator of 
[a/Fe] often used in the literature, of ISM gas (blue, dotted) and 
stars (green, dotted), we find an offset of ~ 0.3 dex with respect 
to the observed [a/Fe]. Note that the size of this offset is depen¬ 
dent on the adopted set of solar abundances. The slope, on the other 
hand, still matches the observed one, implying that the offset can 
be attributed to a constant uncertainty factor in the (massive star) 
yields. 

Motivated by the tight power-law relations shown in Fig. 4, 
we fit the relation between the recycled gas contribution to the SFR 
and ISM metallicity with the following function, including a term 
describing the variation in the relative channel contributions: 


o_ 

Fe 


logi 


SFR let 

SFR 


= O.871og 10 Z gas — 0.40 


+ 0.90, (10) 


where the values of the three free parameters have been obtained 
using least square fitting. Note that the metallicity Z is the average 
mass fraction of metals and is thus independent of the adopted so¬ 
lar value. Similarly, we determine the best-fit relation between the 
recycled gas contribution to the stellar mass and stellar metallicity: 


Jogio 


FI* ,rec 

M* 


0.91 log 10 Z * 


0.28 


+ 0.92. 


(ID 


We show these relations in Fig. 4 for galaxies with M * ~ 10 9,5 
Mg (red, solid line), M * ~ 10 10 ' 5 Mg (blue, dot-dashed line) 
and M, ~ 10 11 ' 5 Mg (purple, dashed line), where we use the 
median values of [a/Fe] gas and [a/Fe]* in stellar mass bins of 
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0.2 dex centred on the respective masses. As expected from Fig. 5, 
the relations at AT* ~ 10 9 ' 5 M 0 and M, ~ 10 10 ' 5 M 0 are con¬ 
sistent, as a result of the median [a/Fe] (of gas and stars) being 
constant for A/* < 10 10 ' 5 M 0 . On the other hand, galaxies with 
AT* ~ 10 11 ' 5 M 0 have SFR lec /SFR and AT* jrec /M* that are 
~ 0.06 and ~ 0.04 dex lower at fixed metallicity due to an en¬ 
hancement in the contribution from massive stars relative to that 
from intermediate-mass stars (reflected by their enhanced [a/Fe] 
abundance ratio). These offsets are somewhat larger than the la- 
scatter in the relation for all galaxy masses (which is set by the scat¬ 
ter at AT* < 10 1() M 0 ), indicating that the variation of the channel 
contributions at AT* > 10 10 ' 5 M 0 significantly impacts upon the 
relation between metallicity and recycling-fuelled star formation 
in high-mass galaxies. It leads to a reduction of SFR rec /SFR and 
AT*, rec /AT» at fixed metallicity that increases with stellar mass, and 
will therefore make any turnover or flattening at the high-mass end 
of the relation between recycled gas contributions and stellar mass 
(as seen in the mass-metallicity relation; see Tremonti et al. 2004; 
Gallazzi et al. 2005; Kewley & Ellison 2008; Andrews & Martini 
2013; Zahid et al. 2014b) more pronounced. We demonstrate the 
useful link that equations (10) and (11) provide between the impor¬ 
tance of gas recycling and observational diagnostics in Section 3.3. 

We note that the parameters of equations (10) and (11) are in¬ 
sensitive to the specific implementation of subgrid processes like 
star formation, stellar feedback and AGN feedback 8 , as for EA¬ 
GLE changing their implementation affects the recycled gas con¬ 
tributions and metallicities in a similar way. Note that this may not 
be true if the metallicity of galactic winds differs significantly from 
the metallicity of the ISM, as might for example happen if metals 
are preferentially ejected (e.g. Mac Low & Ferrara 1999; Creasey 
et al. 2015), or if instead galactic winds are metal-depressed (e.g. 
Zahid et al. 2014a). 


3.3 Dependence on halo and galaxy mass at 2 = 0 

In this section we investigate how the fractional contribution of re¬ 
cycled gas to the present-day SFR and stellar mass of galaxies de¬ 
pends on their halo and stellar mass. Note that, because of the tight 
relation with metallicity that we established in Section 3.2, many 
conclusions that we draw here carry over to the mass-metallicity 
relation. We study both central (Section 3.3.1) and satellite (Sec¬ 
tion 3.3.2) galaxies and, in addition to the total contribution of gas 
recycling, assess the relative significance of the different mass loss 
channels (Section 3.3.3). We also briefly explore how fuelling by 
gas recycling depends on the distance from the galactic centre (Sec¬ 
tion 3.3.4). While we mainly present results from our fiducial Ref- 
L100N1504 simulation, we also show a brief comparison with the 
results from Recal-L025N0752 for central galaxies. 


3.3.1 Gas recycling in central galaxies 

Fig. 6 shows the contribution of recycled stellar mass loss to the 
present-day SFR and stellar mass of central galaxies as a function 
of their mass in the Ref-L100N1504 (red) and Recal-L025N0752 
(purple) simulations. We plot SFRrec/SFR. in the top row and 
AT*, rec /AT* in the bottom row as a function of subhalo mass (left 


8 The adopted IMF is an exception, as it determines the mass and metallic¬ 
ity of gas returned by stellar populations, as well as the relative contribution 
from massive stars with respect to intermediate-mass stars. 


column) and stellar mass (right column). Focusing first on the fidu¬ 
cial Ref-L100N1504 simulation, the general trend in all four pan¬ 
els is that, at masses AT sub < 10 12 ' 2 M 0 or AT* < 10 10 ' 5 M 0 , 
the fraction of the SFR and stellar mass contributed by recycling 
increases with mass. This is the regime where the greater depth 
of the gravitational potential well, as well as the higher pressure 
and density of the ISM and CGM, towards higher masses, make 
it harder for feedback (dominated by star formation) to eject gas 
from the galaxy. As we will show explicitly with a model compari¬ 
son in Section 4, a reduced efficiency of stellar feedback at driving 
galactic outflows enhances the contribution from recycled gas to 
the SFR and stellar mass. This can be understood by considering 
that these winds (if stellar feedback is efficient) are launched from 
the dense star-forming regions with relatively high abundances of 
gas from stellar mass loss. Hence, more efficient winds will prefer¬ 
entially reduce SFR Tec with respect to the total SFR (thereby re¬ 
ducing SFR lec /SFR), whereas in the case of less efficient winds 
this effect will be less (thereby enhancing SFR lec /SFR). 

At the high-mass end, SFR rec /SFR and A/* !re c/AT* turn 
over at AT* ~ 10 10 ' 5 M 0 (AT aub ~ 10 12 ' 2 M 0 ), and then de¬ 
crease and remain constant, respectively, at higher masses. In this 
mass regime, the trend is regulated by the efficiency of the feed¬ 
back from AGN, which becomes stronger in more massive systems. 
Even though this type of feedback is not associated with any replen¬ 
ishment of the ISM gas reservoir (as opposed to feedback from star 
formation, which directly provides the gas for recycling), it does 
have a significant impact on the rate at which galaxies consume 
the enriched ISM gas. If AGN are efficient at launching galactic 
outflows, they preferentially remove or disperse the dense ISM gas 
from the central regions, in which the abundance of stellar ejecta is 
high, thereby reducing SFR rec /SFR and AT*, rec /AT*. 

We note that while AGN feedback regulates the turnover at 
AT* ~ 10 10 ' 5 M 0 , the use of the 30 pkpc 3D aperture also plays 
a role in shaping the behaviour of SFRrec/SFR and AT* !reC /AT* 
at the high-mass end. As we show in Appendix C, SFRrec/SFR 
and AT*, re c/AT* at AT* > 10 10 ' 5 M 0 (A/ aub > 10 12 ' 2 M 0 ) are 
somewhat enhanced, and their slopes become somewhat shallower, 
if an aperture is applied. This is consistent with the fraction of the 
SFR and stellar mass fuelled by recycled gas being larger in the in¬ 
ner parts of galaxies (see Fig. 10). Without an aperture, AT* !re c/AT* 
decreases with halo and galaxy mass (similar to SFR le c / SFR), in¬ 
stead of remaining roughly constant if an aperture is applied. 

At the mass scale of the turnover, the fractional contribution 
of recycled gas to the SFR is at a maximum. A galaxy of this mass, 
AT* ~ 10 10 ' 5 M 0 , is too massive to have effective star formation- 
driven outflows but still too small for AGN feedback to be effec¬ 
tive. Not surprisingly, this mass scale coincides with the peak in 
the galaxy formation efficiency (see fig. 8 of S15), which is con¬ 
sistent with the efficiency of feedback being the main driver of 
SFRrec/SFR and AT* irec /AT». The fiducial EAGLE model indi¬ 
cates that for a Milky Way-like galaxy, which is at the peak of the 
galaxy formation efficiency, 40% of its present-day SFR and 20% 
of its present-day stellar mass is due to the recycling of stellar mass 
loss. 

Because of the tight correlation between recycling-fuelled star 
formation and metallicity, our findings have direct implications for 
the origin of the mass-metallicity relations for ISM gas and stars. 
They imply that the increase in metallicity with stellar mass at 
AT* < 10 10 ' 5 M 0 is due to the decreasing efficiency of stellar feed¬ 
back at driving galactic outflows, while the shape at higher mass is 
governed by the efficiency of AGN feedback (see also Zahid et al. 
2014a; Peeples et al. 2014; Creasey et al. 2015, for discussion on 
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Figure 6. The fractional contribution of gas recycled from stellar mass loss to the SFR (top) and stellar mass (bottom) of central galaxies at z = 0 as a 
function of their subhalo mass (left) and stellar mass (right). We show the results for the fiducial EAGLE model ( Ref-L100N1504: red) and the high-resolution, 
recalibrated model ( Recal-L025N0752\ purple). We only consider subhaloes with a non-zero SFR (top panels) or a non-zero stellar mass (bottom panels). 
The curves show the median value in each logarithmic mass bin of size 0.2 dex, if it contains at least 10 galaxies. The shaded regions mark the 10th to 90th 
percentile ranges. The solid curves become dotted when the subhalo (stellar) mass corresponds to fewer than 100 dark matter (baryonic) particles and become 
dashed (for Recal-L025N0752 only) when there are less than 10 haloes per bin.The contribution of recycled gas to the SFR and stellar mass first increases 
with mass, turns over at M* ~ 10 10 5 Mq (A/ si] |, ~ 10 12 2 Mq), and then decreases or remains constant at higher mass. This trend is regulated by the 
efficiency of the feedback from star formation (AGN) at low (high) masses: galactic winds eject gas from the ISM, where stellar mass loss accumulates, and 
therefore preferentially reduce the SFR and stellar mass contributed by recycling. The black points represent our best estimate of the recycled gas contributions 
to the SFR and stellar mass (for a central galaxy with a Milky Way-like mass: 35% and 20%, respectively), calculated by applying equations (10) and (1 1) to 
the observed gas-phase metallicities from Zahid et al. (2014b) and the observed stellar metallicities from Gallazzi et al. (2005). For M* > 10 111 Mq, these 
estimates agree to better than a factor of ~ 1.6 (0.2 dex) with the median predictions computed directly from EAGLE. 


the relation between feedback and metallicity). Conversely, the dif¬ 
ference between the Ref-L100N1504 and Recal-L025N0752 sim¬ 
ulations in Fig. 6, as well as their (expected) agreement with ob¬ 
servations, should mimic the results for the mass-metallicity rela¬ 
tion (see fig. 13 of S15). Indeed, while Ref-L100N1504 and Recal- 
L025N0752 yield similar trends, they do differ quantitatively by a 
factor of~ 2 (0.3 dex) in SFR lec /SFR and M*, rec /M* at M* ~ 
10 9 Mq (M su b ~ 10 11 Mq). This difference decreases towards 
higher masses, where for M* > 10 9 ' 8 Mq (Af su b > 10 11 ' 6 Mq), 
Ref-LI 00N1504 and Recal-L025N0752 are converged in terms of 
M* <Tec /M* and broadly consistent in terms of SFR roc /SFR (con¬ 
sidering the substantial amount of scatter and relatively poor sam¬ 
pling of the high-mass regime by the Recal-L025N0752 model). 
S15 showed that for M* > 10 9 ' 8 Mq, the metallicities of galaxies 


in Ref-L100N1504 and Recal-L025N0752 agree with the observa¬ 
tions equally well. They agree with the observed gas-phase metal¬ 
licities from Zahid et al. (2014b) to better than 0.1 dex and with 
Tremonti et al. (2004) to better than 0.2 dex, and with the observed 
stellar metallicities from Gallazzi et al. (2005) to within the obser¬ 
vational uncertainties (which are > 0.5 dex at M* < 10 1() Mq 
and smaller at higher masses). For M* < 10 9 ' 8 Mq, on the other 
hand, the metallicities of galaxies in Recal-L025N0752 are in bet¬ 
ter agreement with the observations, from which we conclude that 
the values of SFR rec /SFR and M», re c/A-/* predicted by Recal- 
L025N0752 are more reliable than those predicted by the fiducial 
model. Note, however, that the large systematic uncertainties as¬ 
sociated with the calibration of the diagnostics prevent any strong 
conclusions. In order to limit the number of model curves plot- 
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ted in each figure, from here on we only plot the results from Ref- 
L100N1504 and ask the reader to keep in mind the slight overpre¬ 
diction of SFR lec /SFR and A-f* irec /M, at M* < 10 9 ' 8 M 0 . 

Finally, in contrast to the predictions computed directly from 
EAGLE, which at low masses depend on the adopted numeri¬ 
cal resolution, the relations between gas recycling and metallic - 
ity given in equations (10) and (11) provide a way of estimating 
SFR Tec /SFR and M*, rec /A/*, that is independent of the resolu¬ 
tion. Moreover, these relations are insensitive to the subgrid mod¬ 
els for feedback. We apply the relations to the observed mass- 
metallicity relations from Zahid et al. (2014b) and Gallazzi et al. 
(2005), using the median [O/Fe] from EAGLE in each stellar mass 
bin, to estimate SFRrec/ SFR (triangular points, upper-right panel 
of Fig. 6) and M* jrec /Af* (circular points, lower-right panel of 
Fig. 6) as a function of stellar mass. These estimates agree qualita¬ 
tively with SFR Tec /SFR and A/*, rec /Af* computed directly from 
the fiducial EAGLE model, showing a steep increase with mass for 
M* < 10 10 ' 5 M 0 , followed by turnover and even a slight down¬ 
turn in SFR Tec /SFR at higher masses 9 . Quantitatively, the black 
points are in good agreement with the fiducial EAGLE model for 
M* > 10 10 Mg and with Recal-L025N0752 also at lower masses, 
as expected from the comparison of the mass-metallicity relation 
with the observations presented in S15. If the discrepancy between 
the predicted and observed mass-metallicity relation exceeds the 
systematic error due to calibration uncertainties in the observations, 
then the black points represent our best estimates of the recycled 
gas contributions to the SFR and stellar mass. For a Milky Way¬ 
like galaxy (M* ~ 10 10 ' 5 Mg), we find these contributions to be 
35% and 20%, respectively. 


3.3.2 Gas recycling in satellite galaxies 

Having studied the recycling-fuelled star formation in present-day 
central galaxies, we now compare these with the results for present- 
day satellite galaxies. Fig. 7 shows the SFR and stellar mass con¬ 
tributed by recycling for both central (red; as in Fig. 6) and satellite 
(blue) galaxies, as predicted by the fiducial Ref-L100N1504 simu¬ 
lation. In general, these are broadly similar for centrals and satel¬ 
lites. However, we identify two important differences. Firstly, in 
the left panels, where we show the two ratios as a function of sub¬ 
halo mass, the relations for satellite galaxies are shifted towards 
lower masses relative to those for central galaxies. Satellites lose 
a fraction of their dark matter subhalo mass (but less stellar mass) 
upon infall onto the group dark matter halo as a result of tidal strip¬ 
ping. Hence, this shift illustrates the fact that satellite galaxies live 
in smaller (sub)haloes than central galaxies of similar stellar mass. 
Secondly, in the top-right panel, at a mass scale of M* ~ 10 10 ' 5 
M 0 , satellites show significantly higher SFRrec/SFR (with a me¬ 
dian of ~ 0.5 and a 90th percentile value of ~ 0.85) than cen¬ 
trals, whereas at lower and higher masses this difference is smaller. 
Hence, in the regime where both stellar feedback and AGN feed¬ 
back are relatively inefficient, gas recycling plays a more important 
role in fuelling ongoing star formation in satellite galaxies than in 
central galaxies. For M*, rec /Af*, on the other hand, there is no dif- 


9 Note that, even though the mass-metallicity relation observed by Zahid 
et al. (2014b) does not exhibit a decrease in the metallicity at the high- 
mass end, the recycled gas contribution to the SFR can still show a slight 
downturn, due to the change in the relative contributions from the different 
mass loss channels (as discussed in Section 3.2). 


ference between centrals and satellites, because satellites formed 
the majority of their stars while they were still centrals. 

To get a better understanding of the difference between cen¬ 
trals and satellites, we consider the relation between SFR rec /SFR 
and specific SFR (= SFR/M *, sSFR). Fig. 8 shows this relation 
for centrals (upper panels) and satellites (lower panels) with masses 
10 9 - 5 M 0 < M, < 10 10 ' 5 M 0 (left) and 10 10 ' 5 M 0 < M, < 
10 11 ' 5 M 0 (right), where the histograms at the top compare the dis¬ 
tributions of sSFRs. In order to limit the dynamical range plotted, 
galaxies with SFR/M * < 10~ 12 yr -1 are shown as upper lim¬ 
its. The colour coding indicates the mass of the parent dark matter 
halo, M 200 , in which these galaxies reside. For centrals, M 200 is 
generally closely related to the mass of the subhalo and the stel¬ 
lar mass, whereas for satellites the mass of the host halo they have 
fallen onto is only weakly related to their own mass. For satellites, 
A /200 instead serves as a proxy for the strength of any environ¬ 
mental effects like ram-pressure stripping (Gunn & Gott 1972) or 
strangulation (Larson et al. 1980). 

Focusing first on the centrals in the lower stellar mass bin, we 
see a clear anticorrelation between the fraction of the SFR that is 
fuelled by recycling and the sSFR. This can be explained by the fact 
that sSFR is closely related to the gas fraction (= M gas /(M gaa + 
M t ) with M gaa the ISM mass): higher gas fractions generally cor¬ 
respond to higher sSFRs. Also, since the ISM is comprised of both 
processed and unprocessed gas, whereas stellar mass only provides 
enriched gas for recycling, an enhanced gas fraction typically im¬ 
plies that a greater fraction of the star formation in the ISM is fu¬ 
elled by unprocessed, ‘nonrecycled’ gas (i.e. SFR lec /SFR is low). 
On the other hand, galaxies with a low sSFR and a correspondingly 
low gas fraction will have a higher fraction of their SFR contributed 
by recycling. Considering the centrals in the higher mass bin, we 
see that while part of this relation is still in place, a significant 
fraction lies away from the main relation towards the lower left. 
This is the result of efficient AGN feedback, which suppresses both 
the sSFR (‘quenching’) and the SFRrec/SFR of galaxies, as AGN 
feedback is more important at higher halo masses. This explanation 
is consistent with the enhanced halo masses of the galaxies in this 
regime (M 200 ~ 10 13 M 0 for centrals with SFR/M * ~ 10 -11 ' 5 
yr -1 and SFRrec/SFR ~ 0.2, compared to JV /200 ~ 10 12 M 0 for 
centrals with SFR/M * ~ 10 -10 ' 5 yr -1 and SFRrec/SFR ~ 0.5). 
At even higher central galaxy mass scales than shown in Fig. 8, 
the anticorrelation between sSFR and SFR rec /SFR disappears 
entirely. Instead, the relation transforms into an AGN feedback- 
controlled correlation (although weak due to small number statis¬ 
tics, with a Pearson correlation coefficient of 0.43 for M* > 10 11 0 
M 0 ), such that the galaxies with the lowest sSFRs have the lowest 
recycling-fuelled SFRs. 

Having investigated the mechanism driving the sSFR - 
SFR rec /SFR trends in the two mass regimes for centrals, we now 
consider satellites. In the lower mass bin the anticorrelation be¬ 
tween sSFR and SFRrec/SFR is similar to that for centrals, al¬ 
though the histogram at the top shows that the sSFR distribution for 
satellites has larger scatter towards low sSFRs. This corresponds 
to a population of satellite galaxies at SFR/M » ~ 1CP 11 yr -1 
with fractional contributions of recycled gas to their SFR as high as 
90 — 95%. This high- SFR rec /SFR regime is more frequently pop¬ 
ulated by satellites than similarly massive centrals, reflecting the 
substantially greater scatter in the satellite curves towards high val¬ 
ues of SFRrec/SFR at M* ~ 10 10 ' 5 seen in Fig. 7. As indicated 
by the colour coding in Fig. 8, this population of satellites is hosted 
by relatively massive group dark matter haloes (M 200 ~ 10 14 
M 0 ), implying that their low sSFRs and gas fractions are the re- 
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Figure 7. As Fig. 6, but showing the results for central (red) and satellite (blue) galaxies from the fiducial EAGLE model. The recycled gas contributions to 
the SFR and stellar mass in satellites are broadly consistent with the ones in similarly massive centrals, since the efficiency of stellar and AGN feedback is 
the controlling factor in fuelling star formation with recycled gas. However, in the inefficient feedback regime (M* ~ 10 10,5 Mq), satellites with low gas 
fractions can reach recycling-fuelled SFR fractions as high as ~ 90%, with a median that exceeds the one in similarly massive centrals (see also Fig. 8). 


suit of the cessation of fresh gas infall (either because cooling is 
inefficient or because the satellite’s hot gas reservoir was stripped), 
and/or a (partial) removal of cold gas from the disc (see e.g. Bahe & 
McCarthy 2015; Mistani et al. 2015). Both scenarios lead to a de¬ 
pletion of the ISM gas reservoir and a greater dependence on stellar 
mass loss for replenishing it. 

Finally, we focus on the satellite galaxies in the higher mass 
bin, shown in the bottom-right of Fig. 8. Whereas most similarly 
massive central galaxies that have moved away from the sSFR - 
SFR lec /SFR anticorrelation, have moved towards low sSFR and 
low SFR rec /SFR under the influence of AGN feedback, there 
is still a significant population of satellite galaxies occupying the 
high-SFR lec /SFR region. Inspecting the masses of the BHs resid¬ 
ing in these satellites (not shown) we see that they are significantly 
lower than the masses of BHs in centrals of similar stellar mass. 
We infer that this is again due to the depletion of the satellite ISM 
gas, thereby preventing efficient BH growth. This explains why the 
AGN feedback in these satellites is unable to suppress the recycled 
gas contribution to the SFR. 

We conclude that the SFR and stellar mass contributed by re¬ 
cycling are broadly consistent between central and satellite galax¬ 


ies over a wide range of galaxy masses, because gas recycling is 
governed primarily by the efficiency of stellar and AGN feedback. 
However, in satellites with a stellar mass similar to that of the Milky 
Way, the mass scale at which feedback is least efficient at suppress¬ 
ing star formation, the recycled gas contribution to the SFR often 
exceeds the one in similarly massive centrals (and can even reach 
> 90%), as the depletion of their ISM gas reservoir makes them 
more reliant on stellar mass loss for fuelling ongoing star forma¬ 
tion. 


Our findings are consistent with the observational inference 
that, at a given stellar mass, satellites are more metal-rich than cen¬ 
trals (Pasquali et al. 2012; Peng & Maiolino 2014). We explain 
the origin of their different mass-metallicity relation as a conse¬ 
quence of satellites being subject to environmental processes like 
ram-pressure stripping and strangulation, which prevent the dilu¬ 
tion of the ISM reservoir by metal-poor gas. 
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Figure 8. The SFR fuelled by recycling as a function of sSFR (= SFR/M*), colour-coded by host halo mass (M 200 ), for central galaxies (upper panels) and 
satellite galaxies flower panels) with stellar masses 10 9 5 Mg < M* < 10 10 ' 5 Mg (left) and 10 10 5 Mg < M * < 10 115 Mg (right). The histograms 
at the top show the distributions of the sSFR for centrals and satellites in these two mass bins. Galaxies with SFR/M * < ] 0 12 yr” 1 are shown as upper 
limits. For the centrals, the relation between the recycling-fuelled SFR and the sSFR changes from an anticorrelation at lower mass, which is a result of the 
tight relation with ISM gas fraction, to a (weak) correlation at higher mass, which is driven by AGN feedback. The satellites, on the other hand, show a similar 
behaviour, but retain in both mass ranges a large population of low-sSFR galaxies that rely heavily on stellar mass loss for fuelling ongoing star formation 
(contributing > 90%). 


3.3.3 Contributions from AGB stars, SN Type la and massive 
stars 

To assess the relative significance of the different stellar mass loss 
channels for fuelling star formation in present-day centrals and 
satellites, we show in Fig. 9 the contribution of recycled gas to 
the SFR (top panels) and stellar mass (bottom panels) split into 
the contributions from AGB stars (blue), massive stars (purple) and 
SN Type la (cyan). These are plotted as a function of galaxy stellar 
mass for centrals (left panels) and satellites (right panels) at z = 0. 

From the top panels, we see that AGB stars are of greater im¬ 
portance for fuelling present-day star formation than massive stars 
in all but the most massive central galaxies. Up to 24% (32%) of 
the SFR in centrals (satellites) is fuelled by gas recycled from AGB 
stars, while < 17% (< 20%) is fuelled by gas from massive stars. 
Integrated over cosmic history (as quantified by M*, rec /M, in the 
lower panels), their contributions are approximately equal at all but 
the highest mass scales. This may appear difficult to reconcile with 
the timed mass release from intermediate-mass and massive stars 
for a single SSP presented in Fig. 1, where we showed that mas¬ 
sive stars are the dominant source of (integrated) mass loss for SSP 


ages < 1 Gyr and that massive stars and AGB stars contribute about 
equally at higher ages. Fig. 9 implies that stellar ejecta do not sim¬ 
ply accumulate in the ISM, but that they are removed, either by star 
formation or by outflows, in a way that affects the ejecta from AGB 
stars and massive stars differently 111 . 

The ejecta from massive stars are released almost instanta¬ 
neously compared to those from AGB stars (see Fig. 1 ). This means 
that, before AGB stars start to contribute significantly to the recy¬ 
cled gas, secondary generations of stars will already have formed 
from the massive star ejecta, causing an increasing fraction of these 
ejecta to become locked up in stellar remnants. While they then 
still contribute to the stellar mass of the galaxy, they no longer fuel 
ongoing star formation. Furthermore, a considerable amount of gas 
from massive stars will already have been ejected in galactic winds, 
before the AGB stars shed most of their mass. As a result, the con- 

10 As our findings in Sections 3.3.1 and 3.3.2 already imply, the relation 
between the mass loss from individual SSPs and the contribution of this 
mass loss to the SFR and stellar mass on galactic scales is not straightfor¬ 
ward, since the rate of gas accretion and the efficiencies of stellar and AGN 
feedback depend on the mass of the galaxy. 
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Figure 9. The contribution of gas from AGB stars (blue), massive stars (purple) and SN Type la (cyan) to the SFR (top) and stellar mass (bottom) of galaxies 
at z = 0 as a function of their stellar mass. The results for centrals and satellites are shown in the left and right panels, respectively. The curves and shaded 
regions indicate the medians and 10th to 90th percentile ranges, as in Fig. 7. In general, the gas from AGB stars and massive stars contributes about equally 
to the SFR and stelllar mass in both centrals and satellites. However, there is a slight enhancement in the contribution from AGB stars to the SFR at all but the 
highest mass scales, due to the preferential removal of massive star ejecta by star formation-driven winds and by lock up in stellar remnants. AGB ejecta are 
also responsible for the high SFR Tec /SFR values of some satellites, since these environmentally quenched objects have low sSFRs and cannot accrete gas. 
At the high-mass end, the relative contribution from AGB ejecta declines, because AGN feedback (which is unimportant at low mass) can drive them out even 
in the absence of star formation. 


tribution from massive stars to the present-day SFR is suppressed 
compared to that from AGB stars. 

In general, if stars are older, the surrounding gas will contain 
a higher fraction of ejecta from intermediate-mass stars than from 
massive stars. This means that, at fixed stellar mass, galaxies with 
a low SFR (which includes passive galaxies) are expected to con¬ 
tain enhanced fractions of AGB ejecta, as only newly formed stars 
produce massive star ejecta, while evolved, intermediate-mass stars 
still shed mass in AGB winds. Not surprisingly, we see that the scat¬ 
ter to high values of SFR lec /SFR for satellites with masses 10 10 
Mg < M, < 10 11 Mg (upper-right panel of Fig. 7) is mainly due 
to satellites with large fractions of the AGB ejecta being recycled, 
consistent with their low sSFRs (Fig. 8) 11 . 

11 Despite the shallow, but significant, decrease in the sSFRs with mass of 
galaxies with M* < 10 1(l ° Mg (shown in fig. 11 of S15), we do not see 
an increase in the SFR Iec / SFR from AGB stars relative to SFR Iec / SFR 


Finally, in massive galaxies (M* > 10 10 ' 5 Mg), there is a 
decline of the SFR rec /SFR and M*, re c/Af* contributed by AGB 
stars, while the contributions from massive stars decrease only 
mildly (or flatten). This is consistent with the increase in the [O/Fe] 
abundance ratio, as shown in Fig. 5 of Section 3.2, where we dis¬ 
cussed that this change in the relative significance of the different 
mass loss channels introduces a mass dependence in the relation be¬ 
tween recycling-fuelled star formation and metallicity. We attribute 
this effect to ‘downsizing’, a scenario in which the bulk of the stars 


from massive stars in this mass range. Note that we are now considering the 
ratio of the blue and purple curves in the upper and lower left panels (fo¬ 
cussing on central galaxies). Instead, the relative contribution from massive 
stars, both to the SFR (upper panels) and stellar mass (lower panels), re¬ 
mains approximately constant. This is due to a competing effect, namely 
the decrease in star formation feedback efficiency, which mitigates the pref¬ 
erential expulsion of massive star ejecta. 
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in more massive galaxies has formed earlier and over a shorter pe¬ 
riod of time than in lower-mass counterparts (e.g. Cowie et al. 1996; 
Neistein et al. 2006; Cattaneo et al. 2008; Fontanot et al. 2009). 

The rapid and efficient star formation in the progenitors of 
present-day massive galaxies is however suppressed at later times, 
when these progenitors have grown massive enough for AGN feed¬ 
back to become efficient and the gas cooling rates to drop. More¬ 
over, while winds driven by feedback from star formation will not 
be available to drive out AGB ejecta in quenched galaxies, AGN 
feedback can. Flence, this scenario is consistent with the reduction 
in the contribution from AGB stars, relative to that from massive 
stars, to both the SFR (upper panels) and the stellar mass (lower 
panels) at the highest mass scales shown in Fig. 9. 

We infer that galaxies generally obtain most of their metals 
from massive star ejecta, as these ejecta have 4 — 6 times higher 
metallicity than those from AGB stars, while the fraction of the (to¬ 
tal) ISM mass and stellar mass contributed by massive star ejecta 
is similar to that contributed by AGB ejecta. This holds for metals 
in the gas-phase, but even more so for metals in the stellar phase. 
In both cases, the metal content contributed by massive star ejecta 
increases at the high-mass end (AT* > 10 10 ' 5 Mg) of the mass- 
metallicity relation. This is reflected by the trend of [O/Fe] with 
stellar mass as presented in Fig. 5, and is consistent with the abun¬ 
dance ratio trends observed for early-type galaxies (e.g. Schiavon 
2007; Thomas et al. 2010; Johansson et al. 2012; Conroy et al. 
2014). Our results imply that this a-enhancement of massive galax¬ 
ies is a consequence of AGN feedback. 

3.3.4 Radial dependence of gas recycling 

Having explored the importance of gas recycling on galaxy-wide 
scales, we now briefly investigate how the significance of recycling 
for fuelling star formation depends on the distance from the galac¬ 
tic centre. Note that the robustness of the results that we present in 
this section depends on the ability of the simulation to reproduce 
observed metallicity gradients. In addition, the results are subject 
to numerical uncertainties, like the mixing of metals, which may 
be underestimated by SPH simulations (see Wiersma et al. 2009b). 
This will be investigated in a future paper. In Fig. 10 we plot, sim¬ 
ilar to the left column of Fig. 9, the contribution of gas from AGB 
stars (blue), massive stars (purple) and SN Type la (cyan) to the 
SFR (top) and stellar mass (bottom) of centrals at z = 0, now sep¬ 
arated into gas and stars inside (solid lines) and outside (dashed 
lines) the stellar half-mass radius. This radius, denoted by R 50 , is 
the 3D radius that encloses 50% of the stellar mass bound to the 
subhalo (within the 30 pkpc 3D aperture). It is typically ~ 4 pkpc 
for a AT, ~ 10 1() Mg galaxy. Note that we split both the numer¬ 
ator and the denominator of SFR lec /SFR and M*, rec /M* into 
R < R 50 and R > R 50 . We also plot the total recycled gas con¬ 
tribution (black), which is the sum over the three stellar mass loss 
channels, to the SFR and stellar mass in these two radial regimes. 

Focusing first on the total (black lines), both panels consis¬ 
tently show that gas recycling is more important for fuelling star 
formation (at the present-day and in the past) in the central parts 
of galaxies than in the outskirts. This is consistent with the obser¬ 
vational inference that galaxies grow in an inside-out fashion (e.g. 
Muiioz-Mateos et al. 2007; Patel et al. 2013), with the oldest stars 
residing in the centre and the replenishment of the gas reservoir 
by late infall being primarily significant in the outskirts of the disc 
(owing to its relatively high angular momentum). From the black 
curves in the upper panel we see that, at the peak value, 55 — 60% 
of the SFR inside R 5 q is due to stellar mass loss, compared to only 



M* [M 0 ] 

Figure 10. As the left column of Fig. 9, but split into gas and stars inside the 
stellar half-mass radius R 50 (solid) and outside i?so (dashed). The 10th to 
90th percentile ranges are only shown for inside R 50 . Gas recycling is more 
important for fuelling star formation (at the present-day and in the past) in 
the inner parts of galaxies than in the outskirts. Consistent with inside-out 
growth, the gas in the central regions is comprised of an enhanced fraction 
of AGB ejecta, which is the main driver of the greater contribution of gas 
recycling to the star formation within R 50 . 


35 — 40% outside Rso- As a result, also a higher fraction of the 
stellar mass inside R$o comes from recycling, ~ 30% at the mass 
scale where recycling is most significant, compared to ~ 20 % in 
the outskirts. 

To investigate what drives this radial dependence, we turn to 
the relative significance of the different sources of mass loss. While 
outside R 5Q AGB stars (blue, dashed) and massive stars (purple, 
dashed) contribute about equally to the SFR for all masses, inside 
R 50 the contribution from AGB stars (blue, solid) is significantly 
larger than the contribution from massive stars (purple, solid). As 
discussed in the previous section, the gas around older stellar pop¬ 
ulations contains higher fractions of AGB ejecta. Hence, the dif¬ 
ference between the inner and outer parts reflects the radial age 
gradient of the stars, due to the inside-out growth of the galaxy. 

The drop in the AGB contribution at the high-mass end that 
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we saw in Fig. 9, is present in both the inner and outer parts, but it 
is much stronger near the galactic centre. This is consistent with it 
being due to a lack of intermediate-age stars and the ability of AGN 
to drive out the AGB ejecta even in the absence of star formation. 

The radial variation of SFR rec /SFR and M*, lec /M* is con¬ 
sistent with the negative metallicity gradients observed in local disc 
galaxies (e.g. Zaritsky et al. 1994; Moustakas et al. 2010; Sanchez 
et al. 2014). Our results imply that, although the majority of the 
metals in galaxies (both in the ISM and stars) comes from mas¬ 
sive star ejecta, the relative contribution from AGB ejecta is on 
average larger near the galactic centre than in the outskirts. This 
holds in particular for the ISM metal content of Milky Way-like 
(M, ~ 10 10 ' 5 M 0 ) galaxies. 


4 EXPLORING MODEL VARIATIONS WITH OWLS 

In this section we assess the sensitivity of our results from the EA¬ 
GLE simulation presented in Section 3 to the physical processes 
included in the subgrid model. We do this by comparing a set of 
OWLS simulations (Schaye et al. 2010) in which the subgrid model 
is systematically varied. In particular, we explore variations of the 
feedback from star formation and AGN (Section 4.1), in combina¬ 
tion with metal-line cooling, by explicitly turning on or off a partic¬ 
ular process. The variation of SFR rec /SFR and M*, r ec/AT* with 
galaxy mass in Figs. 6 , 7, 9 and 10 is already suggestive of the im¬ 
portant role played by these feedback processes and we will now 
show this explicitly. We also vary the adopted IMF (Section 4.2), 
as this determines the fraction of the stellar mass that is released 
by a stellar population 1- . The sets of variations are summarized in 
Table 1 and are described in more detail below. 

We start this section with a brief overview of the OWLS 
simulation set-up and the implemented subgrid physics, focusing 
in particular on the differences with respect to EAGLE (for a 
detailed description of the differences, we refer the reader to S15). 
The OWLS simulations were run with a modified version of the 
SPH code GADGET3, but in contrast to EAGLE it uses the entropy 
formulation of SPH implemented by Springel & Hemquist (2002). 
The adopted cosmological parameters, [O m , Ob, Ha, ag,n a ,h\ = 
[0.238,0.0418,0.762,0.74,0.951,0.73], are consistent with 
WMAP 3-year (Spergel et al. 2007) and WMAP 7-year data (Ko¬ 
matsu et al. 2011). The simulations used here were run in periodic 
volumes of size L = 100 /i _1 cMpc. 13 , containing N = 512 3 
dark matter particles with initial mass mdm = 4.1 x 10 s h- 1 M 0 
and an equal number of baryonic particles with initial mass 
mb = 8.7 x 10' ft -1 Mg. The gravitational softening length is 
7.81 h~ x ckpc, limited to a maximum of 2.00 /t -1 pkpc. 

The implementations of radiative cooling and heating, and 
stellar evolution, are nearly the same as in EAGLE. Hence, the 
timed mass release by an SSP, as presented in Section 2.2, is nearly 
identical in OWLS. On the other hand, differences in the OWLS 
subgrid physics include the use of a fixed density threshold in the 
implementation of star formation and the kinetic implementation of 

12 We use the OWLS simulations to perform the model comparison, as this 
suite provides both the extreme feedback variations and IMF variations we 
need. Note that while the stellar and AGN feedback models in OWLS are 
slightly different from those in EAGLE, and that OWLS does not reproduce 
the observed z ~ 0 GSMF, we can still use the OWLS suite to study the 
relative changes in SFR rec /SFR and .VL. rec /A7~ (at least qualitatively). 

13 Note that for OWLS, the box size and particle masses are given in units 
with h,- 1 . 


star formation-driven winds. A fixed fraction f t h of the available 
feedback energy is injected locally, where it directly (by ‘kicking’ 
gas particles surrounding newly-formed star particles) generates 
galactic winds with initial velocity v w and mass loading 77 (fol¬ 
lowing Dalla Vecchia & Schaye 2008). The prescriptions for BH 
growth and AGN feedback (Booth & Schaye 2009), which are only 
included in the OWLS model AGN, employ a Bondi-Hoyle accre¬ 
tion rate with a density-dependent correction term, and a slightly 
lower thermal heating temperature than in EAGLE. 

4.1 Effect of feedback processes and metal-line cooling 

We consider the following set of feedback variations, combined 
with variations in the metal-line cooling, as the latter may also im¬ 
pact upon the efficiency of the feedback: 

• REF is the OWLS fiducial model, which serves as the ref¬ 
erence in the model comparison. It includes radiative cooling and 
heating, star formation, stellar evolution and kinetic energy feed¬ 
back (with fth = 0.40, v w = 600 km s ” 1 and 77 = 2) from 
star formation. Note that it does not include prescriptions for BH 
growth and AGN feedback. All model variations listed below are 
varied with respect to this model. 

• In NOZCOOL the cooling rates are calculated using primor¬ 
dial element abundances, i.e. X = 0.76 and Y = 0.24. 

• NOSN turns off all energy feedback mechanisms associated 
with star formation. However, the mass loss and metal production 
by massive stars, SN Type la and AGB stars are still present. 

• NOSN_NOZCOOL is a combination of the previous two mod¬ 
els. It does not include galactic winds and the cooling rates are 
based on primordial element abundances. 

• AGN includes models for the growth of BHs and feedback 
from AGN. 

Fig. 11 shows the effects of star formation feedback, AGN 
feedback and metal-line cooling on SFR rec /SFR as a function of 
halo mass, M 2 oo 14 - Note that the mass scale corresponding to 100 
dark matter particles, below which the curves are shown as dotted 
lines, is a factor of ~ 58 greater than in the EAGLE fiducial simu¬ 
lation. 

First comparing NOSN (green curve) to REF (black curve), 
we see that the feedback associated with star formation dramat¬ 
ically reduces the contribution of recycled gas to the SFR at all 
mass scales, in particular for M 200 < 10 1- Mg, where the regu¬ 
lation of star formation is largely governed by feedback from star 
formation. Here, the contribution drops from ~ 55% (without star 
formation feedback, NOSN ) to ~ 25% (with star formation feed¬ 
back, REF). The large reduction can be explained if we consider 
the environments from which these star formation-driven winds are 
launched. The outflows originate in the dense ISM, which is the 
environment into which stellar mass loss is deposited. These stellar 
ejecta, which could be recycled into new generations of star for¬ 
mation, are prevented from forming stars if the winds eject the gas 
from the ISM in such a way that it does not return and become 
sufficiently dense again on short timescales. If, on the other hand, 

14 Note that we do not treat central and satellite galaxies separately, but 
instead consider the total amount of star formation taking place within 
the group halo as a whole (without applying an aperture when calculating 
SFi? rec / SFR). Also, while we only present the results for SFR Iec /SFR 
as a function of halo mass, they are consistent with the results for 
A/* r ec /M*, as well as for both ratios as functions of stellar mass. 
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Table 1. Set of OWLS simulations that vary in terms of the feedback implementation (upper section) or the adopted IMF (lower section). From left to right, the 
columns show the model name, whether or not there is energy feedback associated with star formation (SF feedback), metal-line cooling and AGN feedback, 
the adopted IMF, the fraction of kinetic energy available from SN Type II that is used to drive galactic winds (/thX the initial wind velocity (u w ) and the wind 
mass loading parameter (if). 


Name 

SF feedback 

Metal-line cooling 

AGN feedback 

IMF 

fth 

% 

[km/s] 

V 

Feedback variations 

REF 

/ 

/ 

- 

Chabrier 

0.40 

600 

2.0 

NOZCOOL 

/ 

- 

- 

Chabrier 

0.40 

600 

2.0 

NOSN 

- 

/ 

- 

Chabrier 

0.40 

600 

2.0 

NOSN_NOZCOOL 

- 

- 

- 

Chabrier 

0.40 

600 

2.0 

AGN 

/ 

/ 

/ 

Chabrier 

0.40 

600 

2.0 

IMF variations 

1MFSALP 

/ 

/ 

- 

Salpeter 
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a At high pressures ( P/k > 2.0 x 10 6 cm 3 K) the IMF switches from Chabrier (2003) to a top-heavy power-law dN/dM oc M 1 . 
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Figure 11. Comparison of a set of OWLS models to explore the effects of 
star formation feedback, AGN feedback and metal-line cooling on the con¬ 
tribution of recycled gas to the SFR at z = 0. The curves show the median 
(in logarithmic mass bins of size 0.2 dex containing at least 10 haloes) con¬ 
tribution of recycled gas to the SFR as a function of halo mass. The grey 
shaded region shows the 10th to 90th percentile region for the OWLS fidu¬ 
cial model (black curve). The solid curves become dotted when the halo 
mass corresponds to fewer than 100 dark matter particles. We find that 
the efficiency of the feedback associated with star formation (in low-mass 
galaxies) and the feedback of AGN (in high-mass galaxies) determines how 
much gas from stellar mass loss contributes to the SFR. A higher feedback 
efficiency results in a lower contribution from recycled gas. 


star formation feedback is inefficient, then this gas remains in the 
ISM and fuels star formation. The REF model shows an increas¬ 
ing SFR lec /SFR with mass, which becomes increasingly similar 
to the NOSN model. This is consistent with SN feedback becom¬ 
ing less efficient as the depth of the potential well and the density 
and pressure in the ISM increase. Hence, the decreasing efficiency 


of the feedback leads to gas recycling being increasingly important 
for fuelling star formation in more massive systems. 

A notable feature in the curve of the REF model is the sharp 
upturn of SFR lec /SFR at a halo mass of M 200 ~ 10 12 Mq. 
Clearly, the feedback suddenly becomes very inefficient at this 
mass scale. As explained by Dalla Vecchia & Schaye (2012). this is 
due to strong artificial radiative losses in the ISM as the gas (which 
has a high pressure and density in these high-mass systems) gets 
shock-heated by the star formation-driven winds. Kinetic energy 
is thermalised to temperatures at which the cooling time is short 
relative to the sound crossing time, and is quickly radiated away. 
As a consequence, the winds stall in the ISM before they can es¬ 
cape the galaxy, which drives up the value of SFR lec /SFR. For 
a fixed initial velocity of the kinetically implemented winds, this 
effect causes a sharp transition at the mass scale for which artificial 
losses become significant. This results in unrealistic stellar mass 
fractions in haloes of M 200 > 10 12 Mq (Haas et al. 2013) and a 
failure of the model to reproduce the observed GSMF (Crain et al. 
2009). Comparing REF to NOZCOOL (red), we see that the upturn 
becomes less pronounced if the cooling rates are reduced. Lower 
cooling rates, as a result of neglecting metal-line cooling ( NOZ¬ 
COOL), reduce the artificial thermal losses in the ISM and there¬ 
fore enable the feedback to remain more efficient to higher mass 
scales. This results in lower values of SFR Tec /SFR in the most 
massive systems. 

In addition to its impact upon the feedback efficiency, a change 
in the cooling rates also affects the accretion rate onto the galaxy, 
therefore impacting upon SFR rec /SFR in a more direct fashion. 
However, from comparing NOSN (with metal-line cooling; green) 
and NOSN_NOZCOOL (without metal-line cooling; yellow) we see 
that, in the absence of energy feedback associated with star forma¬ 
tion, the effect of changing the cooling rates on SFR Tec /SFR is 
small, especially considering the expected amount of scatter in the 
two relations (from the grey shaded region). Hence, we conclude 
that a change in the cooling rates, as a result of turning metal-line 
cooling on or off, mainly affects the contribution of recycled gas to 
the SFR by changing the (partly numerical) efficiency of the star 
formation feedback implementation. 

Finally, to investigate the effect of AGN feedback on 
SFR rec /SFR as a function of halo mass, we compare the AGN 
model (blue curve), which includes AGN feedback, to the REF 
model. Since BHs live in the dense, central regions of galaxies. 
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where a large fraction of the stellar ejecta are deposited, we again 
expect a low recycling-fuelled SFR if the AGN feedback is strong 
enough to eject gas from the ISM. The AGN model curve indeed 
shows that at masses M 200 > 10 12 Mg, where feedback from 
AGN becomes important, the values of SFR Tec /SFR decrease to¬ 
wards higher masses. At a halo mass of M 200 ~ 10 13 Mg, recycled 
gas contributes only ~ 20% to the SFR. This is in stark contrast to 
the REF model, for which the contribution reaches ~ 65%, indi¬ 
cating the strong impact of AGN feedback on the SFR vec /SFR 
in the regime where the feedback from star formation is inefficient 
and AGN are the main drivers of galactic winds. This highlights 
the importance of including AGN feedback in the subgrid model. 
Qualitatively, we conclude that for massive galaxies the increas¬ 
ing efficiency of AGN feedback towards higher masses leads to gas 
recycling being less important for fuelling present-day star forma¬ 
tion. 


4.2 Effect of changing the stellar initial mass function 

We consider the following variations with respect to the fiducial 
Chabrier (2003) IMF: 

• IMFSALP adopts a Salpeter (1955) IMF, spanning the same 
stellar mass range as the Chabrier IMF used by the fiducial model. 
The corresponding change in the amplitude of the observed Ken- 
nicutt (1998) relation is taken into account. While v w and r/ in the 
implementation of star formation-driven winds are kept the same, 
as is the total wind energy per unit stellar mass (which is propor¬ 
tional to v^rj 2 ), fth is increased to 0 . 66 . 

• DBLIMFCONTSFML14 assumes an IMF that becomes top- 
heavy in high-pressure environments. For P/k > 2.0 x 10° 
cm ' 3 K the IMF switches from a Chabrier IMF to a power-law 
dN/dM oc M -1 . In these environments there is 7.3 times more 
stellar feedback energy available per unit stellar mass to drive 
galactic winds. In this model the additional energy is used to in¬ 
crease the wind mass loading by a factor of 7.3. 

• As the previous model, DBLIMFCONTSFV1618 switches to 
a top-heavy power-law IMF for stars forming in high-pressure re¬ 
gions. However, the additional stellar feedback energy is now used 
to increase the initial velocity of the winds: v w is a factor of x/7.3 
higher than in the reference model. 

Fig. 12 shows the effect of varying the choice of IMF. A com¬ 
parison of IMFSALP (red curve) and REF (black curve) shows that 
adopting a more bottom-heavy IMF like Salpeter reduces the values 
of SFR Tec /SFR over the whole mass range. This was expected, 
because both the total mass and the metal mass released by stel¬ 
lar populations are lower than for a Chabrier IMF (compare Figs. 1 

and Al). In low-mass galaxies (in haloes with masses A /200 < 10 12 

Mg) the contribution of recycled gas to the SFR drops from ~ 25% 
to ~ 15%, while in high-mass galaxies it drops from ~ 65% to 
~ 40%. In the high-mass systems, we expect that the reduction 
is partly due to a reduction of the cooling rates, which causes a 
slight increase of the efficiency of star formation feedback. This 
is the result of the reduced metal mass released by stellar popu¬ 
lations, similar to disabling metal-line cooling in the NOZCOOL 
model (see Section 4.1). However, note that since the total wind 
energy per unit stellar mass is kept fixed when switching from REF 
to the IMFSALP model, the feedback efficiency is affected only by 
a change in the cooling rates. 

On the other hand, the switch to a top-heavy IMF as imple¬ 
mented in DBLIMFCONTSFML14 and DBLIMFCONTSFV1618, 



Figure 12. As Fig. 1 1 , but showing a set of OWLS models with different 
IMFs. We find that adopting a more bottom-heavy (top-heavy) IMF reduces 
(enhances) the contribution from recycled gas to the SFR, but only if the 
feedback energy used to initiate star formation-driven galactic winds is kept 
fixed. If the extra stellar feedback energy from adopting a top-heavy IMF is 
used to either increase the mass loading or the wind velocity, then the recy¬ 
cled gas contributions decrease, showing that other IMF-related effects are 
more than compensated for by the increased efficiency of the star formation 
feedback. 

leads to competing effects. A top-heavy IMF yields more stel¬ 
lar mass loss, but also yields more stellar feedback energy to 
drive galactic winds, which is partly dissipated due the increased 
metal mass loss. Comparing REF with DBLIMFCONTSFML14 
(increased r/) and DBLIMFCONTSFV1618 (increased v w ) enables 
us to determine which effect dominates. We have investigated (but 
do not show here) the effects of increasing the mass loading or the 
initial wind velocity without changing the IMF. We find that if the 
winds efficiently escape the galaxy (in low-mass systems, M 200 
10 12 Mg) the mass loading mainly determines the gas mass that is 
ejected, whereas the initial velocity is of little importance. On the 
other hand, in high-mass systems (M 200 > 10 12 Mg), where the 
artificial radiative losses are high and star formation-driven winds 
are not efficient in escaping the galaxy, we find that increasing the 
mass loading has little effect, as these losses remain too significant. 
Boosting instead the initial wind velocity, alleviates these losses 
and increases the efficiency of the wind, because the wind now ther- 
malises at a higher temperature (Haas et al. 2013). 

This is consistent with the IMF variations, DBLIMFCON- 
TSFML14 (yellow curve) and DBLIMFCONTSFV1618 (purple 
curve), shown in Fig. 12. For DBLIMFCONTSFML14 the values of 
SFR rec /SFR are reduced at low masses and consistent with REF 
at high masses, whereas for DBLIMFCONTSFML14 SFR lec /SFR 
is consistent at low masses and reduced at high masses. Hence, we 
infer that the feedback efficiency is the dominant factor in deter¬ 
mining SFR rec /SFR as a function of halo mass. Despite the fact 
that, naively, we might have expected an increase of SFR lec /SFR 
upon adopting a top-heavy IMF in high-pressure regions, as a re¬ 
sult of the larger fraction of the stellar material available for recy- 
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cling, the decrease of SFR iec /SFR in the DBLIMFCONTSFML14 
and DBLIMFCONTSFV1618 models compared to the REF model 
shows that this effect is more than compensated for by the increased 
efficiency of the stellar feedback. 

From the OWLS model comparisons presented in this section 
we conclude that the efficiency of the feedback from star formation 
and AGN is key for regulating the fuelling of star formation with 
recycled gas. The choice of the IMF sets the total mass of gas that 
is potentially available for recycling, but the feedback efficiency 
determines how much gas recycled from stellar mass loss actually 
contributes to the SFR (and hence the stellar mass) of galaxies. This 
makes the contribution from recycled gas to the SFR and stellar 
mass sensitive to the mass of the galaxy. 


5 SUMMARY AND DISCUSSION 

We have investigated the significance of stellar ejecta as fuel for 
star formation using the Ref-L100N1504 cosmological simulation 
from the EAGLE project. We studied the contribution of gas from 
evolved stellar populations to the SFR and stellar mass, as a cosmic 
average as a function of redshift and within individual galaxies as 
a function of metallicity and galaxy stellar mass at a = 0. We 
treated the galaxies identified as ‘centrals’ separately from those 
identified as ‘satellites’. Since the mass released by AGB stars, SN 
Type la and massive stars was explicitly followed in the simulation, 
we were able to assess the relative significance of these different 
mass loss channels for fuelling star formation. We also explored the 
radial dependence of gas recycling, by comparing the significance 
of recycling-fuelled star formation in the inner and outer parts of 
galaxies. Our results can be summarized as follows: 

• The contribution of recycled gas to the present-day SFR and 
stellar mass of galaxies is strongly, positively correlated with, re¬ 
spectively, the metallicity of the ISM and stars. Therefore, many 
of our conclusions on the role of stellar ejecta in fuelling star for¬ 
mation as a function of galaxy mass and type carry over to the 
mass-metallicity relation. The relations between the contribution 
of stellar mass loss and metallicity do exhibit a slight dependence 
on galaxy stellar mass, as a result of the increasing contribution 
of mass loss from massive stars relative to that from intermediate- 
mass stars to the SFR and stellar mass for M* > 10 10 ' 5 M 0 (Fig. 9, 
Section 3.3.3). We provide the best-fit relations (equations 10 and 
11 ), including a term with the [O/Fe] abundance ratio, which en¬ 
able one to estimate the importance of gas recycling in present-day 
galaxies from the observed metallicity and a-enhancement (Fig. 4, 
Section 3.2). 

• We apply the relations between the SFR contributed by recy¬ 
cling and ISM metallicity and between the stellar mass contributed 
by recycling and stellar metallicity from EAGLE to the observed 
mass-metallicity relations to estimate the recycled gas contribu¬ 
tions as a function of galaxy stellar mass. Since we find these re¬ 
lations to be insensitive to the subgrid models for feedback, ap¬ 
plying them to the observed mass-metallicity relations yields more 
accurate estimates for the contribution of recycling than the di¬ 
rect predictions of EAGLE, provided that the (systematic) uncer¬ 
tainty in the calibration of the observed mass-metallicity relation is 
smaller than the discrepancy between the mass-metallicity relation 
predicted by EAGLE and the observed relation. For central galaxies 
with a stellar mass similar to that of the Milky Way (M* ~ 10 10 ' 5 
M 0 ), which corresponds to the mass scale of the peak in the galaxy 
formation efficiency, 35% of the present-day SFR and 20% of the 
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present-day stellar mass is due to recycled stellar mass loss (Fig. 6 , 
Section 3.3.1). 

• Recycling of stellar mass loss becomes increasingly important 
for fuelling star formation towards lower redshift. At the present- 
day, the fiducial EAGLE model (i.e. as computed directly from the 
simulation) indicates that approximately 35% of the cosmic SFR 
density and 19% of the cosmic stellar mass density is contributed 
by recycling (Figs. 2 and 3, Section 3.1). 

• The fraction of the present-day SFR and stellar mass of cen¬ 
tral galaxies contributed by recycling shows a characteristic trend 
with the mass of the galaxy and its subhalo: for M, < 10 10 ' 5 M 0 
(AT S ub < 10 12 '" M 0 ) the contribution increases with mass, while 
for AT* > 10 10 ' 5 M 0 the contribution turns over and decreases 
with mass (in case of the SFR) or remains approximately constant 
(in case of the stellar mass). We infer that this trend is regulated 
by the efficiency of the feedback associated with star formation (at 
low mass scales) and AGN (at high mass scales). If feedback is 
efficient in driving galactic winds and thereby ejecting gas from 
the ISM, which is the environment into which stellar mass loss is 
deposited, then this will preferentially reduce the SFR and stellar 
mass contributed by recycled gas (Figs. 6 and 7, Section 3.3). 

• The importance of gas recycling for fuelling ongoing star for¬ 
mation in satellite galaxies is broadly consistent with that for cen¬ 
tral galaxies over a wide range of masses, as recycling is mainly 
governed by the efficiency of feedback. However, the fiducial EA¬ 
GLE model indicates that in satellites with a Milky Way-like mass 
the fraction of the SFR contributed by recycled gas significantly 
exceeds the one in similarly massive centrals, and even reaches 
> 90% for satellites with the lowest gas fractions (Fig. 7, Sec¬ 
tion 3.3.2). We infer that this results from a depletion of the ISM 
gas reservoir of the satellite, either due to the cessation of fresh in¬ 
fall or the removal of gas from the disc, which makes them more re¬ 
liant on stellar mass loss for fuelling ongoing star formation (Fig. 8 , 
Section 3.3.2). 

• As a cosmic average, the gas from AGB stars accounts for an 
increasing fraction of the recycled stellar mass loss towards lower 
redshift. At z > 0.4, however, massive stars still provide the ma¬ 
jority of the gas that fuels the cosmic SFR density through recy¬ 
cling. As a result, massive stars dominate the recycling-fuelled stel¬ 
lar mass density at all redshifts. The contribution from SN Type la 
is always small (Figs. 2 and 3, Section 3.1). 

• Within individual galaxies, AGB stars and massive stars con¬ 
tribute approximately equally to the present-day SFR and stellar 
mass of centrals and satellites. The contribution from AGB stars to 
the SFR is slightly enhanced with respect to the massive star con¬ 
tribution at all but the highest mass scales, which results from the 
preferential ejection of massive star ejecta by star formation-driven 
winds and their early lock up in stellar remnants. At the highest 
mass scales (M* > 10 10,5 M 0 ), on the other hand, we find a rela¬ 
tive enhancement in the contributions from massive stars, which we 
attribute to a downsizing effect, with more massive galaxies form¬ 
ing their stars earlier and more rapidly. Their stellar mass therefore 
preferentially consists of massive star ejecta, which are recycled on 
short timescales (Fig. 9, Section 3.3.3). 

• Exploring the radial dependence of gas recycling within cen¬ 
tral and satellite galaxies, we find that recycling is more important 
for fuelling star formation (at the present-day and in the past) in the 
central parts of galaxies (within 77so) than in the outskirts (outside 
R 50 ), which is consistent with the observationally inferred inside- 
out growth of galaxies. We find that the difference between these 
two radial regimes is predominantly driven by the difference in the 
fractional contribution from AGB stars to the SFR (stellar mass). 
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which is significantly higher than (roughly equal to) the one from 
massive stars inside f?so and roughly equal (lower) outside R$o- 
This radial trend directly reflects the negative stellar age gradient 
with increasing distance from the galactic centre (Fig. 10, Sec¬ 
tion 3.3.4). 

Finally, we assessed the sensitivity of our results from the EA¬ 
GLE simulation to the physical processes in the subgrid model 
using a suite of simulations from the OWLS project. The suite is 
comprised of a set of extreme variations of the feedback model, in 
which star formation feedback, AGN feedback and metal-line cool¬ 
ing are switched on or off entirely (Fig. 11, Section 4), as well as 
a set of variations of the adopted IMF (Fig. 12, Section 4). A sys¬ 
tematic comparison of the results shows that while the total fraction 
of the stellar mass that is available for recycling is determined by 
the adopted IMF, the fraction of the stellar mass loss that is actu¬ 
ally used to fuel star formation is controlled by the efficiency of 
the feedback associated with star formation and the feedback from 
AGN, each affecting the galaxy mass regime where the respective 
feedback process regulates the star formation. 

Consistent with previous studies (e.g. Kennicutt et al. 1994; 
Leitner & Kravtsov 2011; Voit & Donahue 2011), our results em¬ 
phasize the importance of modelling the recycling of stellar ejecta 
in simulations of galaxy formation, and the necessity of account¬ 
ing for such gas in assessments of the ‘fuel budget' of present-day 
galaxies. The fractional contribution of recycled ejecta to the SFR 
and stellar mass is not dominant, but it is also not negligible, and 
it extends the gas consumption timescale significantly beyond that 
implied by the ratio of the instantaneous gas mass and star forma¬ 
tion rate of galaxies. 

The relatively small contribution of recycling to the SFR and 
stellar mass of massive galaxies in our simulations is contrary to 
the naive expectation that the establishment of a hot circumgalac- 
tic medium quenches gas infall and renders the galaxy reliant on 
recycling for continued fuelling. Instead, the simulations indicate 
that the ongoing star formation in massive galaxies is sustained 
mostly by unprocessed gas. An interesting route for future stud¬ 
ies will be to explore whether this gas originates in cooling flows, 
or is stripped from infalling satellite galaxies. 
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APPENDIX A: MASS RELEASED BY AN SSP WITH A 
SALPETERIMF 

Fig. Al shows the total (left) and metal (right) mass released by an 
SSP with a Salpeter (1955) IMF (solid lines) and by an SSP with a 
Chabrier (2003) IMF (dotted lines) in the range 0.1 — 100 Mg as 
a function of age for solar metallicity. The colours in both panels 
have the same meaning as in Fig. 1. The (metal) mass loss is lower 
for a Salpeter IMF than for a Chabrier IMF by a factor of ~ 1.5 
over the whole range of SSP ages plotted. The relative contribu¬ 
tions from massive stars and AGB stars are somewhat lower and 
higher, respectively. The metal mass loss from SN Type la is higher 


for a Salpeter IMF, accounting for ~ 8% of the total metal mass re¬ 
leased. In general, adopting a Salpeter IMF instead of a Chabrier 
IMF increases the relative contributions from intermediate-mass 
stars to the (metal) mass loss, as expected for a more bottom-heavy 
IMF. 


APPENDIX B: NUMERICAL CONVERGENCE 

B1 Relation between recycled gas contributions and stellar 

mass 

We test for numerical convergence with respect to resolution of the 
SFR rec /SFR- M* and A/* jrec /M*- M, relations using a set of 
three EAGLE simulations that were run in volumes of size L = 25 
cMpc. We consider both ‘weak’ and ‘strong’ convergence (follow¬ 
ing the nomenclature introduced by S15) by comparing simula¬ 
tions with recalibrated and non-recalibrated subgrid physics, re¬ 
spectively. We use: 

• One simulation with N = 376 s and with the same subgrid 
model parameters as our fiducial L = 100 cMpc, N = 1504 3 sim¬ 
ulation that was used throughout this work. This simulation (de¬ 
noted Ref-L025N376) has the same resolution as the fiducial simu¬ 
lation. 

• One simulation with N = 752 3 and the same subgrid model 
parameters as the fiducial simulation, but with 8 times higher mass 
resolution (Ref-L025N752). 

• One simulation with N = 752 3 and with a recalibrated set of 
subgrid model parameters for star formation feedback, AGN feed¬ 
back and the accretion onto BHs, in order to improve the match 
with the observed 2 ~ 0 GSMF at this 8 times higher mass res¬ 
olution ( Recal-L025N752 ). In short, the recalibration corresponds 
to a change in the density dependence of the stellar feedback effi¬ 
ciency parameter ft h, such that the feedback efficiency is increased 
in higher-density gas while keeping the average f t h roughly equal 
to 1. This is done in order to compensate for the increase in cooling 
losses, which arise as a result of the locally higher gas densities that 
are resolved in the higher-resolution model. 

A comparison of these three EAGLE simulations is shown in 
Fig. Bl, where we plot SFR lec /SFR (top) and M*, re c/Af* (bot¬ 
tom) as a function of stellar mass. We show these relations for both 
central galaxies (left) and satellite galaxies (right). Since our con¬ 
clusions about resolution convergence are broadly consistent be¬ 
tween centrals and satellites (although with somewhat poorer sam¬ 
pling for the latter), we will focus the discussion below on central 
galaxies. 

Comparing the fiducial resolution simulation ( Ref 
L025N0376 ; red) to the two higher-resolution simulations, 
with ( Recal-L025N0752\ purple) and without (Ref-L025N0752; 
blue) recalibrated subgrid feedback parameters, we infer that 
for central galaxies with M* < 10 9 ' 8 Mg, SFR lec /SFR and 
Af*, rec /M* as a function of stellar mass are not numerically 
converged at the fiducial resolution. The ‘strong’ convergence 
is somewhat better than the ‘weak’ convergence. At M* ~ 10 9 
Mg, SFRrec/SFR and M, |rec /Af* in Ref-L025N0376 are almost 
a factor of 2 (0.3 dex on a logarithmic scale) higher than in 
Recal-L025N0752. This is not surprising considering the level of 
agreement between Ref-L100N1504 and Recal-L025N0752 for the 
mass-metallicity relations, where the latter is in better agreement 
with the observations (see fig 13. of S15). 

At masses M* > 10 9 ' 8 Mg, the relation between M*, reC /M* 
and stellar mass is fully converged (both ‘weakly’ and ‘strongly’). 


MNRAS 000. 1-23 (2015) 


22 M. C. Segers el al. 




Age [yr] Age [yr] 


Figure Al. The cumulative fraction of the initial mass (total: left panel; in the form of metals: right panel) that is released by an SSP as a function of its age, 
adopting a Salpeter (1955) IMF (solid lines) or a Chabrier (2003) IMF (dotted lines). The curves show the contributions from AGB stars (blue), massive stars 
(purple) and SN Type la (cyan), as well as the total (metal) mass ejected by the SSP (black) for solar stellar metallicity. The total (metal) mass loss is lower 
for a Salpeter IMF than for a Chabrier IMF by a factor of ~ 1.5. Adopting a Salpeter IMF increases the relative contributions from AGB stars and SN Type 
la, which have intermediate-mass progenitor stars. 


while SFR rec / SFR as a function of stellar mass shows substan¬ 
tial overlap between Ref-L025N0376, Ref-L025N0752 and Recal- 
L025N0752. Due to the small box size, however, SFR Tec /SFR. 
is not well-sampled around AT* ~ lO 10 ' 5 Mg, the mass scale at 
which SFR lec: /SFR reaches a maximum in our fiducial L = 100 
cMpc model (see Fig. 6). 


B2 Relation between recycled gas contributions and 
metallicity 

Fig. B2 shows the numerical convergence test of the SFR lec /SFR- 
Z gas (left) and AT„, rec /AT*- Z , (bottom) relations for central galax¬ 
ies, comparing the fiducial EAGLE model (Ref-LIOON1504\ red) 
and the high-resolution, recalibrated model ( Recal-L025N0752\ 
purple). While Recal-L025N0752 spans a metallicity range that 
is shifted towards somewhat lower values with respect to Ref- 
L100N1504, due to its smaller box size and 8 times higher mass 
resolution (we select stellar masses corresponding to at least 100 
gas particles at each resolution), the SFR rec /SFR- Z gas and 
AT*, rec /AT*- Z * relations are converged with resolution over the 
whole metallicity range probed here. Where Recal-L025N0752 and 
Ref-L100N1504 overlap, their medians agree to better than 0.05 
dex in SFR vec /SFR and to better than 0.04 dex in AT*, rec /AT*. 

One might wonder whether the secondary dependence on a- 
enhancement (hence, implicitly on stellar mass; see Fig. 5) af¬ 
fects the convergence of these relations. Flowever, the dependence 
on stellar mass only becomes significant for AT* > 10 10 ' 5 Mg, 
which is the regime where Ref-L100N1504 and Recal-L025N0752 
are converged (in terms of Af*, r ec/AT* and ZF) or at least broadly 
consistent (in terms of SFR lec /SFR and Z gaa ). As a consistency 
check, we repeat the calculation of SFR lec /SFR and AT*, re c/AT* 
by applying the relations between recycling and metallicity to the 


observed mass-metallicity relations (as done in Section 3.3.1) at 
higher resolution using Recal-L025N0752. We find agreement with 
the results from Ref-L100N1504 to better than a factor of ~ 1.07 
(0.03 dex) over the whole stellar mass range. 


APPENDIX C: EFFECT OF USING A 3D APERTURE 

Fig. Cl shows the effect of using a 30 pkpc 3D aperture on the con¬ 
tribution of recycled gas to the SFR (top) and stellar mass (bottom) 
in central galaxies at 2 = 0 as a function of their subhalo mass 
(left) and stellar mass (right). We compare results from the EAGLE 
fiducial model, Ref-LIOON1504, with (red) an without (purple) the 
aperture. Recall that the aperture only applies to galaxy properties, 
hence the subhalo mass is not affected. 

While the majority of the star formation takes places within 
the central 30 pkpc, causing the effect of the aperture on the total 
SFR to be small, there is still an enhancement in SFR lec /SFR (up¬ 
per panels) inside the aperture compared to its value over the whole 
galaxy. This enhancement is consistent with recycling-fuelled star 
formation being more important in the central parts of galaxies (see 
Fig. 10). The effect is significant for AT* > 10 11 Mg (AT su b > 
10 12 ' 5 Mg) and increases with mass up to a difference of a factor 
of ~ 1.5 (~ 0.18 dex) in the left panel, whereas in the right panel 
the effect is smaller (up to a factor of ~ 1.2, or ~ 0.08 dex) due 
to the simultaneous decrease in the stellar mass if an aperture is 
applied. 

Similarly, in the lower panels, AT,, re c/AT* is enhanced if an 
aperture is applied. In this case, the effect of using an aperture is 
that the decrease of AT*, re c/AT* with mass for AT* > 10 10 ' 5 Mg 
(AFsub > 10 12 ' 2 Mg) becomes a flattening at a roughly constant 
value (~ 22% instead of 17% at AT* ~ 10 11 ' 5 Mg). 
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Figure Bl. Test for numerical convergence of the SFR rec /SFR- M* (top) and M*, rec /M*- M* (bottom) relations (presented in Figs. 6 and 7) for central 
(left) and satellite (right) galaxies at z = 0, comparing an EAGLE simulation with the fiducial resolution {Ref-L025N376\ red) to two EAGLE simulations 
with 8 times higher mass resolution, with {Recal-L025N0752\ purple) and without ( Ref-L025N0752\ blue) recalibrated subgrid feedback parameters. All three 
simulations were run in volumes of size L = 25 cMpc. The curves and shaded regions have the same meaning as in Fig. 6, but for clarity we only show the 
10th to 90th percentile range for the Ref-L025N376 and Recal-L025N752 simulations. For galaxies with M* > 10 9,8 Mq, SFR rec /SFR and M*, rec /M* 
as a function of stellar mass are reasonably well converged. For galaxies with M* < 10 9 ' 8 Mq, the fiducial EAGLE simulation likely overpredicts the SFR 
and stellar mass contributed by recycling, by at most a factor of ~ 2 (~ 0.3 dex) at M* ~ 10 9 M©. 
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Figure B2. Test for numerical convergence of the SFR rec /SFR- Z gas (left) and M* )rec /M*- Z * (right) relations (presented in Fig. 4) for central galaxies 
at z = 0, comparing the fiducial EAGLE model ( Ref-L100N1504\ red) and the high-resolution, recalibrated model {Recal-L025N0752\ purple). We only 
consider galaxies with stellar masses corresponding to at least 100 gas particles, at the respective resolution. In the left panel we only consider subhaloes with 
a non-zero SFR. The curves show the median value in each logarithmic metallicity bin of size 0.05 dex, if it contains at least 10 galaxies. The shaded regions 
mark the 10th to 90th percentile ranges. The SFR rec /SFR- Z gas and M*, rec /M*- Z* relations are converged at the fiducial resolution over the whole 
metallicity range. 
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Figure Cl. The effect of using a 30 pkpc 3D aperture, comparing results from the EAGLE fiducial model, Ref-L100N1504, with (red) and without (purple) 
applying the aperture in calculating galaxy properties. The curves show the fraction of the SFR (top) and stellar mass (bottom) contributed by recycling for 
central galaxies at z = 0 as a function of their subhalo mass (left) and stellar mass (right). The curves and shaded regions have the same meaning as in Fig. 6. 
The effect of using an aperture is a change in the slope of the two recycled gas fractions at M* > 10 11 Mq (M su b > 10 12 ° Mq), which becomes somewhat 
shallower. For the fraction of the stellar mass contributed by recycling this results in a roughly flat trend instead of a decrease with subhalo and stellar mass, 
hence mitigating the effect of increasing AGN feedback efficiency. 
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